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Chapter 1 


INTRODUCTION 

Rising fuel costs have encouraged weight reductions of commercial 
transport aircraft. Weight reductions of approximately 20-30 percent 
may occur with the replacement of existing aluminum construction with 
advanced composite materials. Composite materials offer both high 
strength- and sti f fness-to-weight ratios. Also, the directionality of 
these composite materials allows designers to tailor the properties of a 
structure and, thus, design more efficiently. Composite materials have 
additional advantages over metallic materials in their tension-tension 
fatigue and corrosion properties. 

Although composites offer many inherent advantages over metals, 
there are still many problems which must be solved before they can be 
used for primary load carrying structure. Past research has uncovered 
shortcomings of composites which are not found in metals such as delam- 
ination, free edge effects, and poor impact damage tolerance. Much work 
has been done on the response and failure of composite materials sub- 
jected to in-plane loadings. However, little experimental work has been 
done on composite materials under combined membrane and bending loads. 

A combined stress state occurs in the skin of a fuselage which is inter- 
nally pressurized. This Internal pressure could be carried most effi- 
ciently in a uniform membrane state of stress if the skin were free to 
expand uniformly in the radial direction. Internal stiffeners necessary 
to support mechanical loads applied to a fuselage restrain the skin 
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locally and prevent uniform radial expansion consistent with a membrane 
state. The pressure-i duced deformation of the skin in a typical bay of 
a stiffened fuselage is illustrated in Figure 1 and is often referred to 
as "pillowing." 

The product of the large pressure-generated membrane tension and 
change in slope of the skin's reference surface, which results from non- 
uniform radial expansion, produces a geometrically nonlinear pressure 
response. This nonlinearity is included in this study and is essential 
for accurate response prediction. The degree of radial constraint also 
introduces bending strains which occur in the skin adjacent to the 
stiffeners. The bending strains, which are a maximum aC the panel's 
edge, combine with the membrane tension strains to initiate failure in 
the skin. 

The primary structural configurations of composite pressure vessel 
research are thick unstiffened shells which carry very high pressures. 
Much literature has been published on this research area due to the 
development of high strength fibers and filament winding techniques. 
Analysis of these structures using geometrically linear membrane theory 
is sufficiently accurate for design. Thin stiffened shells such as a 
fuselage require a different procedure. Bending in addition to membrane 
action must be included in the analysis of these structures. Refer- 
ences 1-7 report on studies of fuselage-type structures under Internal 
pressure. Portions of the material presented in this thesis are docu- 


mented in Reference 8 
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Figure 1.- Test specimen’s relationship to an internally pressurized fuselage. 
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Reference 1 reports on a one-dimensional nonlinear elastic membrane 
analysis for long flat or slightly curved metal panels. This work 
illustrates the nonlinear nature of the problem but does not address the 
severe bending gradients which occur in the fuselage skin at the stiff- 
ener support. 

Flugge (Ref. 2) addresses a number of problems which arise in the 
internal pressurization of an airplane fuselage. The importance of 
the nonlinear response is illustrated by a membrane example similar to 
that described in Reference 1. Local bending stresses induced by the 
restraint of the stiffeners on the fuselage skin are not predicted with 
the membrane approximation of Flugge. Neither can these local bending 
stresses be correctly examined with a linear axis ymme trie approximation 
of the pressurized fuselage as used by Flugge. In this approximation 
the skin and longerons are lumped together to obtain an equivalent 
"smeared" orthotropic shell with the individual rings modeled dis- 
cretely. Skin bending adjacent to a ring will not be properly predicted 
with this approach and bending adjacent to a longeron is not allowed by 
the analysis. 

Additional axisymmetric analyses are includea in References 3 
and 4. Williams (Ref. 4) included in his analysis the important geo- 
metric nonlinear term which accounts for the product of the in-plane 
axial load and slope in the radial equilibrium equation. However, the 
importance of the nonlinear effect of this term is lost when "smeared" 
longeron theory is used which makes the skin effectively very thick and 
the contributions of the nonlinear term small. This approach therefore 
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does not accurately address the "pillowing” which results in local skin 
bending. 

Wang (Ref. 5) modeled the fuselage as a shell connected to discrete 
internal rings and stringers. The analysis was linear; and it was 
assumed that the interaction between elements consisted only of normal 
stresses, thus excluding shear. Series solutions were used to satisfy 
the equilibrium equations for each individual component. Displacement 
compatibility between the individual elements couples these elements and 
allows for the determination of the interaction forces. This analysis 
for isotropic materials does address "pillowing" effects but geomet- 
rically linear behavior is assumed. Linear behavior would become more 
appropriate as the skin thickness increases or as the stiffnesses of the 
frames and longerons decrease. 

Formulas are given in Reference 6 for stresses in an isotropic skin 
adjacent to individual frames and longerons. Although the details of 
the derivation are not discussed nor a reference given, nonlinear terms 
are apparent in the formulas given. 

Reference 7 reported axial and circumferential strain distributions 
for a typical bay of a composite stiffened circular cylindrical shell 
under internal pressure. These results were generated from a Lockheed 
in-house computer program. The analysis developed by Wang (Reference 5) 
has been extended to include composite skins. The analysis is again 
linear with numerical results given for a fuselage skin laminate at an 
ultimate ground test condition of 17.63 psi. 
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References 1-7 all deal with pressurization of a fuselage but none 
show how the geometric nonlinearities effect the severity of the bending 
gradients in a composite skin adjacent to a stiffener. Only Reference 6 
has experimental pressurization results. This work is on large scale 
stiffened adhesively bonded aluminum fuselage structures and does not as 
such focus on the response of the skin. For panels with longitudinal 
cracks in the center of the bay, the crack propagated until it met a 
frame and turned and ran parallel to the frame along its intersection 
with the skin. For aluminum pressurized panels, the high stresses which 
occurs in the local bending gradient are relieved by material yielding. 
For composite panels, strain relief must come primarily from other 
sources such as transverse matrix cracking, local fiber failure, and 
delamination. 

In the edge bending gradient region, interlaminar stresses become 
large and must be examined because of inherent transverse strength 
weaknesses of composites. The majority of references on interlaminar 
stresses is directed towards the free edge problem of a composite 
laminate. References 9 and 10 examine the distribution of interlaminar 
stresses near the supported edges of composite circular cylinders under 
internal pressure. Linear shell analysis was used or. a layer-by-layer 
basis to determine the axial distribution of interlaminar stresses for 
clamped or simply supported axial boundary conditions. Because linear 
theory was used in the above studies, the same response characteristics 
would be found at any pressure. The response character will vary with 
pressure when geometrically nonlinear behavior is considered. 



7 


A one-dimensional approach is used in the present study to concen- 
trate on the nonlinearities and local skin bending important in a pres- 
surized fuselage skin. The skin is approximated by an infinitely long 
shallow cylindrical panel. This geometry approximates typical fuselage 
skin bays which are long in the axial direction relative to their cir- 
cumferential arc length. In addition to the in-plane circumferential 
stresses, interlaminar stress distributions are examined. Nonlinear 
equations of elasticity, which include geometric shallowness and 
Donnell's approximations, are integrated to determine the interlaminar 
stresses. Past research has not examined these stresses with geometric 
nonlinearities considered. 

Experiments were conducted to verify the analysis and identify the 
ultimate failure pressures and modes. Comparison between analysis and 
experiment points out the range of applicability of the analysis. The 
test article is a shallow cylindrical panel of the same approximate 
dimensions as a typical fuselage skin panel. Under internal pressure 
the radial deflections of the fuselage skin are symmetrically dis- 
tributed along lines perpendicular to the frames and longerons as shown 
in Figure 1. Because of this symmetry the slope of the deflection curve 
normal to the boundary is zero. To approximate this slope condition, 
curved graphite-epoxy and aluminum panels were tested in a fixtuv? with 
clamped edges. The outward radial deflections at the stiffeners of an 
actual fuselage were not duplicated in the experiment. The panels were 
constructed of 4, 5, 8, or 16 plies of unidirectional graphite-epoxy 
tape to illustrate a range of responses and failures. 



Chapter 2 


EXPERIMENTAL PROGRAM 


2.1 Test Specimens 

The materials used in this study include 2024-T3 aluminum and 
commercially available 0.005-inch-thick unidirectional Thornel 300 
graphite fiber tapes preimpregnated with 350°F cure Narmco 5208 thermo- 
setting epoxy resin referred to as T300-5208 graphite-epoxy. The 
aluminum panels were cut to size from flat sheets which had nominal 
thicknesses of (J.020 and 0.040 inches. Unidirectional preimpregnanted 
tape was layed up on a smooth cylindrical surface with a 60 inch radius 
at specified orientations to form uncured laminates. These panels were 
cured in an autoclave following the resin manufacturer's recommendations 
and then inspected with ultrasonic C-scan. Typical lamina properties 
are 19.0 msi fcr the longitudinal Young's modulus 1.89 msi for the 

transverse Young's modulus ms * f° r t * ie in'Plane and shearing 

deformation moduli and G i 3 > 0*60 msi for the shearing deformation 

modulus g 23 » and 0.38 for the major Poisson's ratio v^* The panels 
and their stacking sequence are shewn in Table 1. The aluminum 
specimens are designated specimens Al and A2 and the graphite-epoxy 
specimens are designated G1 through Gil. 

After trimming the panels to an approximate size of 23 inches by 
11 inches, bolt holes were drilled around the panel's edge for securing 
the panel to the fixture as shown in Figure 2. A special drill fix- 
ture was used for drilling thirty-two 0.25-inch holes into each panel. 
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Table 1 . Test Specimens 



Number 

Specimen of Plies 

Thickness, 

in. 

Laminate Stacking 
Sequence 


A I 

1 

1 

U.UZUJ 

2024-T3 Aluminum 

A2 

1 

0.0389 

2024-T3 Aluminum 

G1 

4 

0.0206 

[,45] s 

G2 

4 

0.0207 

h45] s 

G3 

5 

0.0251 

[ ±45/9$] s 

G4 

5 

0.0247 

[ ±45/515] s 

G5 

8 

0.0387 

[90/+45/0] 

s 

G6 

8 

0.0397 

[90/0/+45] 

s 

G7 

8 

0.0400 

[90/C/±45] 

G8 

8 

0.0912 

^2s 

G9 

16 

0.0901 

[±45/±45/90 2 /0 2 ] 

G10 

16 

0.0880 

r90 2 /±45/±45/0 2 ] 

Gil 

16 

0.0893 

[90 2 /0 2 /±45/±45] 
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The drill fixture consists of a top and bottom part. The trimmed panels 
were placed between the two curved surfaces of the drill fixture for 
drilling the holes. Force was applied to the top and bottom parts of 
the drill fixture with C-clamps to force the panel Into an untwisted 
60 inch radius configuration and to prevent any movement during the 
drilling operation. 

The drill fixture design had an important effect on the circumfer- 
ential distance between bolt holes of different thickness panels. The 
drill fixture design provided for constant opening angles between the 
arc-wise position of the bolt holes drilled in all the panels. However, 
it did not compensate for the various panel thicknesses. Thus differ- 
ences occur in the middle surface arc lengths between the bolt holes 
used in clamping the traight edge. The middle surface arc length dif- 
ference between the thickest panel (0.080 in. nominal) and the thinnest 
panel (0.020 in. nominal) was predicted to be approximately 0.004 in. 
Although this difference may seem like a minor variation in the middle 
surface arc length, it will be shown later that small changes in the arc 
length for shallow panels can result in very large changes in the middle 
surface radius. 

Fiberglass tabs of uniform thickness (approximately 0.1 in.) were 
used to reinforce the bolt holes against bearing failure along the 
straight edges of the panels. These flat tabs measured 1.5 in. in width 
and 20 in. in length and were drilled separately and bonded onto the 
convex surface of the predrilled composite panels. 
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To measure the bending gradient away from the clamped edges, the 
panels were in°trumented with strain gages. Back-to-back strain gages 
were bonded to the panels from the edges to the center along the circum- 
ferential and axial center lines of the panel. Because of the rapid 
strain variation at the panel's edge, strain gages with short gage 
lengths were located as close as possible to the clamped boundary to 
determine the best pointwise estimate of the edge strains. Gages with 
an effective gage length of 0.015 Inches were used close to the edge. 
Since the strains were almost constant outside the edge bending boundary 
layer, larger gages could be used away from the edge. These strain 
gages had an effective gage length of 0.187 inches. 

2.2 Test Apparatus 

The test apparatus consisted of the test fixture, pressure source, 
and instrumentation. A cross sectional drawing of the assembled test 
fixture is shown in Figure 3. The test fixture consists of .three compo- 
nents: test frame, strain gage lead feed-through panel, and top clamp- 
ing bars. The test frame was machined from a solid steel block and has 
a 60-inch radius machined surface around the edges where the curved 
panels were clamped. The strain gage lead feed-through panel allows 
recording a maximum of 16 strain gages bonded to the panel concave sur- 
face. The feed-through panel allows the strain gage signals to be 
transmitted out of the pressurized interior of the test fixture. The 
feed-through panel also has an inlet port for the pressurizing media, a 
bleed valve, and a pressure transducer port. Six bars were used for 
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Figure 3.- Cross section of test specimen in test fixture. 
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clamping the specimen to the fixture. These include two straight bars 
and two curved bars, rectangular in cross section, which were used to 
clamp the panel to the ter,t f.&me. Thirty-two 0.25-inch bolts secured 
these bars and the test panel to the test frame. These bolts were 
torqued to 100 in-lb. Two larger bars, L-shaped in cross section, 
allowed application of very large normal forces to the straight edges of 
the panel in an attempt to minimize in-plane slippage of the panels in 
the clamped support. The bolts which passed through the flange of the 
larger bars are tightened down on the smaller bar as shown in Figure 3. 

A photograph of a failed test specimen assembled in the test fixture is 
shown in Figure 4. The L-shaped bars are not shown in the photograph. 

During the course of the experimental program various methods were 
used to pressurize the curved panels. Pressure sources Include the use 
of low pressure shop air (to 50 psi), city water (to 60 psi), an 
1800 psi bottled nitrogen source, and a hydrostatic tester. A pressure 
regi lator between the pressure source and the interior of the test panel 
was used to control the pressure. 

After the manufacturing and drilling of the composite curved pan- 
els, their shapes deviated from the desired 60-inch radiu6 right cir- 
cular cylinder. Therefore, geometric imperfection equipment was used to 
measure the initial shape of the curved panels. A DCDT (direct-current 
displacement transducer) was fixed to a trolley which rolled on 
straight, flat rails. Scans of the panels were made in the two direc- 
tions parallel to the supports to determine the true geometry of the 


panels. 
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Figure 4.- Test specimen clamped in test fixture. 
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The test procedure consisted of the following steps. A panel was 
placed in the test fixture and the clamping bars were torqued symmetri- 
cally. Next, the panels were surveyed with the geometric imperfection 
measuring equipment. The imperfection measuring equipment was removed 
and DCDT's were positioned normal to the panels surface at selected 
points along the axial and circumferential centerlines to monitor dis- 
placements of the panels to the applied pressure. Pressure was 
increased until failure of the panels occurred. Digital strain, 
displacement, and pressure data were recorded on a data acquisition 
system. 



Chapter 3 


ANALYSIS 

In this chapter the energy Method is used to derive two sets of 
geometrically nonlinear equilibrium equations for the response of 
cylindrical panels to internal pressurization. For a set of nonlinear 
strain-displacement equations, elasticity equations are derived in 
Section 3.1. Using the same set of nonlinear strain-displacement equa- 
tions and explicit assumptions for the displacement variation as a 
function of the thickness coordinate z, shell equations are derived in 
Section 3.2 which account for through-the-thickness shearing deforma- 
tions. In Section 3.3 the two-dimensional shell equations are reduced 
to a one-dimensional set of ordinary differential equations in the 
circumferential coordinate by assuming that the stress resultants are 
independent of the axial coordinate. The reduction to a one-xlimensiona 1 
set of equations is appropriate for the response of a panel which is 
long in the axial direction. Closed form solutions to the nonlinear 
ordinary differential equations are obtained which include transverse 
shear deformation effects and twist-curvature coupling of laminated 
composite panels. The one-dimensional solution provides details of the 
bending gradient response near the straight edge of the panel. In 
Section 3.4 through-the-thickness shear and normal stresses are obtained 
in this bending boundary layer that contains the bending gradient. 

These stresses are determined by substituting the in-plane stresses and 
radial deflection from the one-dimensional shell solution into the 
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elasticity equations derived in Section 3.1. The stresses in the 
elasticity equations are also assumed independent of the axial 
coordinate to be consistent with the one-dimensional shell solution. 
Finally, in Section 3.5 a strength of materials solution is presented 
for the nonlinear response of an internally pressurized cylindrical 
membrane which is long in the axial direction. 

3.1 Three-Dimensional Elasticity Equations 

Elasticity equilibrium equations will be used to determine the 
transverse stresses in Section 3.4 after the shell solution is obtained 
in Section 3.3. However, these equations are developed before the shell 
equations since an elasticity formulation has fewer assumptions, and 
hence, is less restrictive than a shell formulation. 

Figure 5 shows the middle surface or reference surface of the 
circular cylindrical panel. The shell coordinates shown are the axial 
coordinate x, circumferential coordinate 9, and the thickness 
coordinate z. The origin is at the center of the middle surface such 
that -L < X < L, -a < 0 < a, and -t/2 < z < t/2, where 2L is the 
length of the panel, a is the semi-opening angle, and t is the 
thickness. The radius of the circular arc on the middle surface is 
designated by a. 

Let u, v, and w designate displacements in the axial, circum- 
ferential, and thickness directions, respectively, of a material point 


in the panel. The normal 

strains 

are denoted by 

e xx* 

e 86’ and e zz" 

and the engineering shear 

strains 

are denoted by 

Y 0z* 

I'xz- “ nd >x0- 
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The strains are assumed to be small, and the rotations are assumed to be 
moderately small. On a surface parallel to the middle surface, the 
rotations of the axial and circumferential line elements out of the 
tangent plane are assumed to be larger than the rotations about the 
normal, such that rotations about the normal are neglected. Since the 
panel is shallow, the contribution of the circumferential displacement 
to the rotation of the circumferential line element out of the tangent 
plane is neglected (Donnell-Mushtari-Vlasov approximation). Sanders 
(Ref. 11) presents kinematic relations for the surface strains in 
general curvilinear coordinates under these assumptions. Specialized to 
cylindrical coordinates Sanders kinematic relations are 

du 1 ' 5w 
£ xx 2\dx/ 


J. 'f- l/i- dw ^ 
£ ee = r ae + r + 2Vr 56' 


(3.1) 


Sv i 5u 5w 1 &w 
^x0 bx + r 50 + 5x r 30’ 

where r = a + z is the radius of a surface parallel to the middle 
surface. The transverse kinematic relations are assumed to be given by 
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&v 1 5w 
^9z dz r &0 

&w [ du 
^y.z &x &z 


(3.2) 


_ _Qw 
£ z z 5z * 

The panel is thin, such that a » t > |z|. Thus in equations (3.1) and 
(3.2) the radius r is replaced by the radius a of the middle 
surface. 

The potential energy for the panel is 

V = U - W, (3.3) 


where U is the strain energy and W is the work of the prescribed 

surface tractions. Body forces are neglected. Let a , o OQ , and 

xx yy 

0 denote the normal stress components in the panel, and let x_ , 

Z Z 02 

and x denote the shear stress components. The strain energy 
for a lineal elastic thin shell is 



-t/2 -a /-L 

J_ t /2 j- a J-i/ £xx ° xx + C °® a " + £zz 


o 

zz 


9z 0z 


xz 


x + Y „ x ^ dx ad9 dz 
xz 'x0 X0 J 


(3.4) 


The prescribed surface traction components at x = ±L are denoted 

a , x and x . The overbar notation is also used for the pre- 
XX X0 xz 

scribed surface traction components at 0 = ±a and z = ±t/2. The pre 
scribed tractions are treated as dead loads in the energy ormulation. 
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Although hydrostatic pressure is not a dead load, it is conservative, 
and may be approximated by a dead load potential (Ref. 12) for a shallow 
panel. Thus, the work done by the prescribed tractions is 



(3.5) 


where the internal hydrostatic pressure p is equal to ~ a zz at 
z =* t/2. 

Equilibrium equations and boundary conditions for the thin shallow 
cylindrical panel are obtained from stationarity of the potential 
energy, i.e. , 6V - 0. Combining equations (3.3), (3.4), and (3.5) with 
the approximate kinematic relations (3.1) and (3.2), the first variation 
in the potential energy, after integrating by parts, is 



3 





a 

zz 


z=-t/2 


( 3 . 6 ) 
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Since the variations 6u, 6v, and 6w are completely arbitrary, equa- 

tion (3.6) can vanish only if the coefficients of each variation vanish 
individually. From this reasoning the three equations of elasticity 
for equilibrium in the x, 8, and z direction are 


bo , 14 

xx _1_ 

5x a 


r t bx 

XU . xz 


bx 


x8 J_ 


bx 


bx 


bx 


+ — 


+ — 


ae + 

bz 


5t - 

88 

8z 

be 

bz 

bx. 

bo 

8z 

zz 

be 

bz 

b / °86 

l bw\ 

be \ a 

56/ 


(3.7) 


, „ . bw 

- - + -r- 0 . 


a 88 bx \ xx bx 




Tne nonlinearities appear in the z-direction equilibrium equation only, 
which is t.ie third equation in (3.7). 

The vanishing of the first variation also leads to boundary condi- 
tions on the surfaces of the panel. Since the equilibrium equations are 
used to estimate the transverse stress components, only the conditions 
on the upper and lower surfaces are needed. Thus, conditions on the 
surfaces at x = +L and 8 ■ + a are omitted. At z • ±t/2 the 
following are prescribed 


‘‘xz or u ’ 

x or v, (3.8) 


or w 
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3.2 Two-Dimensional Shell Equations 
3.2.1 Kinematic Relations 

The assumptions of geometric shallowness, small strains, and moder- 
ate rotations were used in the derivation of the elasticity equations. 

In this section a shell theory will be derived from the variational 
principle (3.6) by assuming the z-direction dependence of the displace- 
ments. Normals to the reference surface before deformation are assumed 
to remain straight and unchanged in length but not necessarily normal to 
the deformed reference surface. This assumption allows threugh-the- 
thickness shear deformation of the Mindlin type (Ref. 14) to occur which 
is important for composite materials where the ratio of shearing defor- 
mation modulus to in-plane extensional modulus is low (1/20-1/30). This 
assumption implies that the in-plane displacements u and v vary 
linearly through the thickness. Let u° and v° represent the dis- 
placements u and v on the middle surface (z * 0), and let and 

y designate the rotations of the middle surface normal about the 
9-axis and x-axis, respectively. Then the displacement field for the 
shell theory has the form 

u(x,6,z) * u°(x, 9) + z T x (x, 9) 

v(x,9,z) “ v°(x, 9) + z fglx, 9) (3.9) 


w(x, 9,z) * w°(x, 9) 
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Substituting the displacements (3.9) into the strain-displacement 
relations (3.1) and (3.2) results in the following expressions 


c 

e = e + z T 

XX XX XX 


r * r + 7 r 

00 ee 00 


Y 


X 0 




Y 


9z 



(3.10) 



e 

zz 


o 

e 

zz 


0. 


in which the middle surface strain-displacement relations are 


2 

o = 5u° _1_ ( 5w° i 

E xx 5x 2 V 5x 



— dv° + iL. _1_ / J_ 5w° 

a 96 a 2 l a 50 


2 


(3. 1 1) 


o _ 5v° J_ 5u° ftw C J_ 5w° 

Y x0 dx a d9 dx a 60 ' 
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and where Che transverse strains and rot Cion gradients are 


o 



z 


i aw 

v + — -==- 

e a ae 




r 

XX 



• '!Is 

96 a ae 


' x 0 


e + i x 

ax a as 


(3.12) 


3.2.2 Shell Equilibrium Equations 

The assumptions for the displacements (3.9) are substituted into 
the first variation of the potential energy (3.6) and explicit integra- 
tion in the thickness coordinate z from -t/2 to +t/2 is carried 
out. In this process weighted integrals of the stresses in z occur. 
These integrals are interpreted as stress resultants and stress couples 
in shell theory. The stress resultants and couples, defined per unit 
arc length on the reference surface, are 
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After performing the integration through the thickness, and using the 
stress resultant definitions (3.13), the vanishing of the first varia 
tion of the potential energy leads to the equilibrium equations 
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3N . ar: 0 

x J_ x 9 

&x a 3x 


I SN n 3N 

JL 0 + x_e 

a 30 3x 


~ ( N ) 
3x x 3x / 


+ 


dx 


+ 


* lT xz (n/2> - V t/2)l - 

* IV+t/2) - l fe (-t/2)J - 

. l5 + i5_fle *) _,*« 

a a 30 a 30 / + a 30 

ill. j*L \ d_ f ^x0 3w 
a 30 ^ x0 3x / 3x \ a 30 


0 


0 


(3.14) 


+ lo zz (+t/2) - o zz (-t/2)] - 0 

»x 1 » X 0 t 

36 4 a 30 ‘ Q x + 2 It xz (t/:) + \ z { ~ t/2)] " 0 

i*e »x0 r tf 

a 30 dx ' Q 0 + 2 lT 0z (t/2) + T 0 z ( ‘ t/2)] = °- 

The shell boundary conditions, which result from the integrated form of 
the 1 l rst variation of the potential energy, consist of five conditions 
along each edge. At x « ±L the following are prescribed 
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v V o -0 

N * N or u = u , 
x x 


M ■ M or ¥ = ¥ , 

XX XX 


N X0 = *X0 or v ° = *°’ 


M X0 = M X0 or *6 = V 


^ N’ ^ 

n , M 5w x9 - o-o 

Q X - N X to + T Te ‘ Q X or “ 


(3.15) 


at 9 * ±a the following are prescribed 


M V ° “° 

N xe ' \e ° r u ‘ u ’ 


f, x9 ’ M xe or *« ■ V 


— 0 — O 

N q = N g or v = v , 


M 0 = M 0 or ’ 0 - V 


x N 0 3w - o -O 

r si *x8 h ” 5 e or w * k ' 


(3.16) 


The shell equilibrium equations (3.14) can be obtained from the 
elasticity equilibrium equations (3.7) by integrating the latter equa- 
tions through the thickness. For example, if the last ■'uation in 
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equations (3.7) is first Integrated in z from — t / 2 to t/2, and 

secondly it is multiplied by z and integrated from -t/2 to t/2, 

then the first and fourth equations in equations (3.14) are obtained 

when the resultant definition in equations (3.13) a r e used. This fact 

is significant for the estimations of the through-the-thickness stress 

components x , x_ , and a from equations (3.7) when stress com- 

ponents o , o A , Tq, , and deflection w are assumed to be given by 
r x Q tfx 

the shell theory. 

3.2.3 Shell Constitutive Equations 

The panels in this study are fabricated by laminating plies of a 
unidirectional fiber-reinforced composite material. It is assumed each 
ply is homogeneous, linear elastic, and orthotropic with respect to the 
fiber (1), transverse (2), and through-the-thickness (31 directions. 

The fiber direction with respect to the positive x-axis in each lamina 
is given by the angle <J> shown in Figure 5. Thus, in the shell coordi- 
nates x, 9. and z, the material appears to be monoclinic with a plane 
of symmetry perpendicular to the z-axis. In addition it is assumed that 
the through-the-thickness stress a ? is negligible with respect to the 
in-planc stresses a^, o^, and x^q, such that o z is set to zero in 

the constitutive equations. Since the strain e ? also vanishes (see 
equations (3.10)), the constitutive equation for the z-direction normal 
stress and strain is neglected. For the in-plane stresses and strains, 
the lamina constitutive equations are 
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r 'n 
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/ 0- 

V _ 

6 

f 

CD 

X 
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V J 



Q 11 Q 12 Q 16 


Q 22 Q 26 


symn;. 


<66 


f. 

I X 

S c e ) ’ 


(3.17) 


vH 


in which the elements Q ..., Q 66 are the transformed reduced 

stiffnesses determined by the elastic moduli Ej, E 2 , v i 2 * an< ^ ^ 12 ’ 
and the angle $. (See Ref. 13, pp. 46-51.) The lamina constitutive 
equations for the through-the-thickness shear stresses and strains are 


X (tz 


xz 
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— 



r “N 

C44 

C 45 
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V 
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"45 

'55 


Y xz 



__ 


L J 


(3.18) 


where 


'44 


cos" <t + sin 0 G^ 


= sin $ cos <t> (G 2 ^ - G^) 


(3.19) 


C,^ = sin 0 G 01 + cos $ G n , 


23 


13’ 


and Gj and G 23 are the shear moduli in the 1-3 and 2-3 planes, 
respectively. 

The shell constitutive equations are obtained by substituting the 
strains (3.10) into equations (3.17) and (3.18), and then substituting 
these results into the definitions for the stress resultants (3.13). 
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Since the panels tested are balanced and symmetric laminates, the shell 
constitutive equations for this important class of laminates are 


A 11 A 12 


A 12 A 22 


0 0 


0 0 


0 0 0 D, 


0 0 0 D, 




. O 

< Y ::e 


k l A 44 0 


0 k 2 A 55 


where the laminate extensional stiffness (Ajj) and the laminate bending 
stiffness (D^) are determined from the individual lamina stiffness and 
the stacking configuration 


fh/; 

. ,D . . ) = 

J 1J J-h /: 


Q. (l,z )dz (i,j = 1,2,6) 


A. , = 




(l.j ' 4,5) 
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In equation (3.21) kj and k ? are shear correction factors, which are 
introduced to account for the nonuniformity of the through-the-thi ckness 
shearing stress distribution. 

F jrmulas for the shear correction factors are given by Whitney 

(Refs. 15 and 16) for laminated anisotropic plates under static loading. 

For homogeneous isotropic plates various values of the shear correction 

2 

factors have been used (e.g. , 5/6, u / 1 2 , and 2/3; see Ref. 17). The 
shear correction factor serves to change the effective value of the 
laminate shearing deformation modulus. It is assumed that the shear 
correction factors are unity in this analysis, since the values of 
shearing deformation moduli are difficult to measure and are not well 
known. 

3.3 One-Dimensional Shell Equations 

The solution of the two-dimensional shell equations is very diffi- 
cult and requires the use of a numerical computer code. In transport 

aircraft design the axial lengths of fuselage panels are usually greater 
than the circumferential arc lengths or widths. A simpler analysis is 
possible for the central section in long panels if it is assumed that 
the static response to internal pressure is independent of x. Thus 
stresses and strains depend only on 9. The solution to the shell 
equations based on this assumption is developed in this section. 

As reported in Reference 1, W. D. Douglas in 1916 argued that the 
end effects associated with the curved-end boundaries for internally 
pressurized fabric panels are confined to approximately a half panel 
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width in the axial direction from each curved edge. On this basis an 
analysis which is independent of the axial coordinate x is applicable 
for | x | < L - aa. The axial length of the central region where a one- 
dimensional analysis is applicable is examined in more detail in 
Section 4. 1 . 


3.3.1 Shell Equations 

Assuming the stress resultants and couples are independent of the 
x-coordinats the shell equilibrium equations (3.14) simplify to 


1 

a 



+ 



(t/2) 


t ("t/2)] 
xz 


0 


a 7T * 'V t/2) - - 0 


N + + 1 dQ e 

‘x.2 a a 50 \a 59 / ad0 
5x 


1 5 5w \ N ’x9 5 1 2 w 

a 59 ' x 0 5x I a 6x59 


(3.22) 


+ [ o (t/2) - o (-t/2) ] - 0 
zz zz 


dM 




XZ 


1 dM 9 r 

— - Q G + i lT 0z (t/2) - T 9z (-t/2) 1 - 0 


a d 9 



3 b 


3.3.2 General Displacement Field 

The assumption that the stress resultants and couples are indepen- 
dent of the x-coordinate implies the strains and rotation gradients are 
independent of the x-coordinate as well, since the two sets of variables 
are linearly related by the constitutive equations (3.20) and (3.21). 

For this general case, however, the displacements and rotations are not 
independent of the x-coordinate. To determine the explicit dependence 
of the displacements and rotations on x, the strain-displacement 
equations (3.11) and (3.12) are integrated with respect to x under the 
restriction that the strains are functions of 0 only. Integrating 
equations (3.12) gives 

w°(x,9) * f^O) “ ~2 C jX ♦ C 9 x0 + C^x 

Y (x, 9) = f . ( 0) + C x (3.23) 

X 4 1 

T q (x, 0) = f 5 ( 0) - C 2 x, 

in which C j , C 2 , and are arbitrary constants, and f-j, f^, and 

^ are arbitrary functions of 0. Using the functional form of 
w°(x, 0) in the first of equations (3.23) in the process of integrating 
equations (3.11) gives 
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o 12 3 1 • 9 

u°(x,e) * fj(e) + c 1 f 3 (e) x - -g- CjX + y c 1 c 2 x“e 

- j aCj0 2 x + i C j c 3 * 2 + (C 4 - C 2 C 3 ) x0 + fc 5 - j C 2 j x 

(3. 24) 

v°(x,0) = f 2 (e) - - c 2 f 3 (0) x + ^ c 1 c 2 x 3 

1 / 1 2\ 2 1 21 2 
+ 7 i C l - 7 C 2 ' x e - 7 C 2 x9 - 77 V ' c 3 xe + V’ 

where C 4 , C^, and are arbitrary constants, and and f 2 are 

arbitrary functions of 0. 

Using linear elasticity equations Lekhnitsk.il (Ref. 18) develops 
general displacement equations for an anisotropic cylindrical body in 
which stresses are independent of the axial coordinate x. Equivalent 
u° and v° displacements are linear and quadratic in the x-coordinate , 
respectively. Lekhnitskii 's equivalent w° displacement is quadratic 
in the x-coordinate. The difference between the displacements derived 
here and Lekhnitskii 's results is due to the use of nonlinear strain- 
displacement relations and the assumption of a linear displacement vari- 
ation in the z-coordinate. 

3.3.3 Application of Special Boundary Conditions 

Special boundary conditions were applied to the general displace- 
ment relations (3.23) and (3.24). These conditions were selected 
because they are needed to model accurately the experimental boundary 
conditions. They are applied along the circumferential edges of the 
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panel (9 * a and 0 ■ -a) and place further restrictions on the dis- 
placement functional form. The boundary conditions are 


u°(x,±a) 


EX 


v°(x,±a) 


+v 


w°(x,±o) «• 0 (3.25) 

^(x.ta) = 0 

T 0 (x,±a) = 0 

where e and v are prescribed data. Application of these boundary 
conditions reduces the general displacements (3.23) and ( 3 . 2 4) to 


u° (x, 9) = f j ( 0) + ex 
v°(x, 0) = f 2 ( 6) 

w°(x, 9) = f 3 ( 0) (3.26) 

1 (x.0) = f.(0) 
x 4 

¥ e (x,e) = f 5 ( e) , 
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where 


f j (±a) = 0 
f 0 (±a) = +v 

* 0 

f . (±a) = 0 
u 


(3.27) 


f^(±a) = 0. 

3. 3. A Reduction of the Field Equations 

The one-dinensional shell equations are simplified further for the 
particular problem of Interest. In addition to the boundary conditions 
(3.25), the applied loading in the experiment is an internal pressure 
with no surface shear tractions applied to the panel. Hence, for this 
loading the equilibrium equations (3.22) become 
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(3.28) 
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Considering the functional form of the displacements in equa- 
tions (3.26) the strain-displacement relations in equations (3.11) and 
(3.12) reduce to 


o 



^x0 a d0 



_1_ dw° 
a d 0 



r = 0 

xx 
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00 


' X 0 
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Us 

d 0 

df 
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d 0 


(3.29) 


From the first equilibrium equation in equations (3.28) N^q is 
spatially constant. The constitutive equation for N^g in equa- 
tions (3.20) then implies the middle surface shear strain y° is 

x 9 

spatially constant. For a spatially constant middle surface shear 
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strain, the third equation in (3.29) combined with the boundary condi- 
tions on function fj(6) in equations (3.27) requires fj(0) to vanish 

for all Q. r hus N = v° - 0 for all 0. Equations (3.29), with 

x0 'x0 

the addition of fj equal to zero, imply the constitutive 
equations (3.20) reduce to 


N X * A I1* ♦ A i2 c e 


N a ■ A 12 E + A 22 c 6 


«x ■ D 12 r 09 + D 16 r xe 


■ D 22 r ee + D 26 r xe 


‘'xe ' D 26 r 09 + D 66 r x9' 


The constitutive equations (3.21) for the shearing deformation resul- 
tants remain unchanged, and are repeated below for convenience with 

2 2 

shear correction factors kj and k ? assumed to be unity. 
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Using the equilibrium equations (3.28), with N * 0 for all 6, 


In combination with the strain-displacement equations (3.29) and 



constitutive equations (3.30) and (3.31), the mathematical formulation 
reduces to 
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Equations (3.32) to (..36) constitute an eighth order system of ordinary 
differential equations for the dependent variables Nq, w°, T , 
and v°, with p and e as prescribed loads. They are to be solved 
subject to the eight boundary conditions 


v° (ta) 


♦v 


(3.37) 


w°( + a) * 0 


(3.38) 


¥ 0 (±a) - 0 


(3.39) 



* 3 


t Ua) - 0, (3.40) 

x 

which are obtained from equations (3.26) and (3.27) where v is a 
prescribed displacement. 

3.3.5 Method of Solution 

Equations (3.32) to (3.36), subject to boundary conditions (3.37) 
to (3.40), are a nonlinear two-point boundary value problem. The cir- 
cumferential stress resultant Ng is spatially constant as the 
equilibrium equation (3.32) indicates. Since Ng is spatially 

constant, equations (3.33) to (3.35) appear to be linear in w°, f , 

0 

and t . However, Ng depends nonlinearly on w° as shown in equa- 
tion (3.3b). This structure of the boundary value problem permits the 
following solution procedure. Using the methods of solution for linear 
ordinary differential equations with constant coefficients, equat- 
ions (3.33) to (3.35) are solved for w°(0), f (0), and t (0), subject 

6 x 

to boundary conditions (3.38) to (3.40). These solutions will depend on 
the "coefficient" Ng. The solution for w°(0) is then substituted 
into equation (3.36) and indefinite Integration on 0 is perfromecl 
using the fact that Ng is spatially constant. The unknown constant of 
integration and the unknown value of Ng are determined by boundary 
conditions on v° in equation (3.37). 

By neglecting the bending terms (these have coefficients ^22* 

Dgg, and ^26^ equations (3.32) through (3.36) obtained apply to the 
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corresponding membrane shell. Subject to boundary conditions (3.38), 
the membrane solution is 


w°(0) 

P 


a 

2 



(3.41) 


¥ 


0p 


( 0 ) 



(3.42) 


¥ ( 9 ) = 0 . 

xp 


(3.43) 


Notice that in the nonlinear membrane theory it is possible to satisfy 
boundary conditions on the deflection w°(9). In linear membrane shell 
theory the boundary conditions on the deflection cannot be satisfied in 
general. The fact that nonlinear membrane theory can satisfy physically 
reasonable boundary conditions whereas linr ir membrane theory cannot, 
was pointed out by Bromberg and Stoker (Ref. 19). 

A homogeneous solution of equations (3.33) through (3.35) has the 

form 




(3.44) 
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which, when substituted into equations (3.33) through (3.35), leads to 



where \ - — . Nontrivial solutions of equation (3.45) require the 

a 

determinant of coefficients to vanish, which results in the character- 
istic equation 
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The roots of the characteristic equation are 0, 0, ±Xj, and ±X 0 , 

where 



,22 

b + c i 


/ 2 o 4 

V(b l - c ) 


2 2 2 
+ 4 e b c 


!(1 - 


e 2 ) 


(3.48) 


The four nonzero roots given by equation (3.48) are associated with 
bending boundary layers at 6 = ±a, which decay exponentially from the 
edges towards the center of the panel. Approximate expressions for the 
roots may be obtained by performing an asymptotic expansion in the small 
parameter e which is less than unity for the materials examined. 

These approximate expressions reveal more about the nature of the roots 
and are given as follows 




= b‘ 


b 2 + 


b 2 ^ 


0(e 4 ) 


-2 

x 2 



4 

+ 0(e ) 


(3.49) 


The parameter e is zero for specially orthotrcpic laminates (D 26 * 0) 

and the roots X^ and X 2 uncouple. For small e the root X^ is 

associated with the rotation Y (circumferential shear and bending), 

0 

and the root X 0 is associated with the rotation Y (axial shear and 
2 x 

the twisting moment). Fo, the laminates and pressures considered here 
the root X 0 is an order of magnitude larger than the root Xj. Thus, 
the twist-curvature bending boundary layer effects decay at a much 
faster rate from the clamped edges than the circumferential bending 


boundary layer effects 
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The complete homogeneous solution is 



where Cj, Cn, . .., c^ are arbitrary constants, and the eigenvectors 
are the solution of 



The principle of superposition does not hold for equations (3.32) 
to (3.36) because the equations are nonlinear. However the total solu- 
tion for w°, y , and 4* may be obtained by adding the membrane 
8 x 

solution (3.41-3.43) and the homogeneous solution (3.50) where the con- 
stant Ng which occurs in both solutions has not been specified yet. 
The solution is exact once and v° are determined from equations 

(3.36) and (3.37). The six constants for the total solution are deter- 
mined from the six boundary conditions given by equations (3.38) to 
(3.40). Writing the exponential functions in terms of hyperbolic 
functions, the solution subject to the boundary condition is 
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f° ( 0) 


pa - N. 


/ -2 -2 




\ S. /\-2 -2 

■ s - h] 



a..a 55 ] 


D 22 A 55 + D 66 A 44 J_ 
A 4 4 A 5 5 X 


j \ cosh - cosh X^a 

-3 ) sinh X, a (3.52) 

V 1 


pa ~ N e \ I \ \ \ ( °26 D 22 D 66 - 

v n 9 / Vi ~ ^ 2 / A ^ A ^5 


> cosh X2 6 - 

cosh X2C1 

h ' N eiAr 

1 sinh 

X ? a 

A » e A 2 



V 0) 




sinh X.j9 
sinh X^a 


! 253 _ 

\ D 66 V 5lnh X 2 9 

-2 sinh X 9 a 
z L 


+ 




V 0 ) 



D 26 

sinh 

CD 

sinh 

X 9 0 

A 55 

sinh 

\ a 

sinh 

V 


(3.53) 


(3.54) 


The dependent variables w°, T , and ¥ were not evaluated 

D X 

numerically using the functional forms shown in equations (3.52) to 
(3.54). For large values of Xj and the hyperbolic functions 

exceeded the maximum magnitudes permitted on CDC 170 series machines at 
NASA Langley. Since the hyperbolic functions occur in the numerator and 
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denominator of the terras In equations (3.52) to (3.54), it is possible 
to rewrite the devision of two large numbers in a more numerically 
efficient manner. To accomplish this, the hyperbolic functions are 
written in exponential form and then manipulated to make the denominator 
approximately equal to one. A typical example of this procedure is 
shown below. 


, . 0 A , 8 ~A . 6 

cosh A, 9 1 1 

1 e + e 

■■ — . r - 

sinh AjO AjOl - AjO 
e - e 


-AjO 


-Ajd 

e 


Aj(e-a) “Aj (8+a) 

e + e 

-2a , a 
1 - e 


cosh Aj9 Aj(8~a) -Aj(8+a) 

■ e + e 

sinh Ajd 


(3.55) 


Near the edge 9 = a the second term in equation (3.55) is negligible 
with respect to the first such that 


cosh Aj0 Aj(0-a) 

■ - ■ ■ - — • 0 

sinh AjO 


(3.56) 


and near the edge 9 = -a the first term is negligible with respect to 
the second term such that 


cosh A i 9 -Aj(6+a) 

■ ■ - - m 0 

sinh A j a 


(3.57) 


Let $ * 0 + a so that at the edge 9 * -a, $ = 0, and <J> increases 
toward the center of the shell. In terms of the coordinate $, equa- 


tions (3.52) to ( 3 . 54 ) are approximated by 
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w°( $) 


—2 —2 
pa - N Q \ / X> ^ 


\ 2 ' ^ 


P ?6 D 22 D b6 - + P 2 2 A 5 5 + D 66 A 44 j_ 

A 44 A 55 1 A 44 A 55 l 




(-\* , 
e - 1 


pa - N 


-2 -2 
X 1 X 2 
-2 -2 / 
X 1 “ 


D 26 " D 22 D 6b - 


A 44 A 55 


S 


D„ „A C c + D . .A , , 
22 55 66 44 J_ 

A 44 a 55 X. 




-x, <t> 
e - 1 




0 


(3.58) 


V*> 


pa - N 


I? % \ 


N “ 


66 

^55 


^5 . -2 

, °66 2 /’V 

- e 


*4 


- \ 2 

>°66 V-V 

“ * 

X 1 


/pa - n a ' 

+ ! — — 1 ( (J> — ct ) 


\ N 


(3.59) 


e 


? (♦) 

X 


x 2 t2 
e W A i 


pa - N \ I ^ \ D 


26 


“ X*> 4 1 X, 


* 55 L 


e - e 


(3.60) 


The circumferential resultant N 0 is an unknown constant at this 
point. The value of it and the circumferential displacement v°(0) are 
determined by substituting w°(0) given by equation (3.52) into equa- 
tion (3.36), integrating on 0, and applying the two boundary conditions 



51 


piven bv equation (3.37). The result of this procedure for Ng as a 
function of the internal pressure parameter R is the equation 


A(N q ) R 2 + R + C(N n ) = 0 


in which 


Pa - N, 


K i ( i - 

A 1 \ tanh X, a . , 2 

1 1 sinh \.a / 

\ 1 / 


a \ 


4 X„ tanh \ a 2 

^ 2 2 sinh 


*1 K 2 ( X 1 _ *2 \ 
2 _ 2 1 tanh X a tanh X a 

*1 ~ \ \ 1 2 / 


K. / X. a 
1 / 1 




Xj \tanh Xj o J X^ \ tanh X^a 


- 1 6 




Xj 1 tanh X t a X^ j X 0 I tanh X^a X^ I 3 


l N e ' A i2 e ) a v 
a 22 
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and 


I \ \ V U 26 ~ D 22 D 66 -2 A 55 D 22 * A 44°66 

-2 \ A . ,A rr X 1 A, , A , 




44“55 


44 55 


(3.66) 


' t2 \ 

i N *2 ' 

\ [26 " D 22°66 -2 

A 55 D 22 * 

V.i io 

1 -2 -2 

\ A , ,A C c \ 

A , , A » , 

-2 la 

\ X 1 - v 

44 55 

44 55 

X 2 / (3.67) 

equation 

(3.61) are readily 

obtained by 

assuming a value 


of Nq and solving for R, and hence p, as the roots of the resulting 
quadratic equation. In the cases examined in this study the pressure 
associated with one root of equation (3.61) was positive and the 
pressure associated with the second root was negative. The negative 
pressure solution to equation (3.61) was disregarded since it is not 
consistent with an internally pressurized shell. These solutions 
relating the pressure p to the circumferential stress resultant Ng 
are then used to determine the other unknowns of this nonlinear shell 
response problem. 

The exponential terms in equations (3.58) to (.3.60) decay away from 
the boundary ( <t> = 0) because of their negative exponents. A boundary 
decay length may be defined as the distance from the boundary for an 
exponential term to decay to a small value e. After selecting a value 
for t a boundary decay length (£) may be determined 


- Xa <t> 
e = e 
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or 


l 



(3.68) 


where 



A typical value for e might be 0.01 and two values of X are given by 
equation (3.48). Therefore, two boundary decay lengths are evident in 
the dependent variables given by equations (3.58) to (3.60). One decay 
length is associated with twist-curvature effects ( Xj ) and the other 
with circumferential bending effects (X 9 ). 

3.3.6 Klrchhof f-Love Approximation 

In this section additional kinematic restrictions will be placed on 
the panel displacement fields. Normals to the undeformed surface will 
be assumed to remain normal to the deformed surface. This assumption 
means the shearing deformation strains y° and y° vanish, which 
imply the rotations in equation (3.12) become 

* - - (3. u9) 

x dx 

* - - 1 

9 a d0 


(3.70) 
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Thus the displacement fields in equation (3.9) have the following form 
for the one-dimensional analysis 

u(x,6,z)-ex (3.71) 

v(*i,8,z) - V°(e) - z - 4r-(e) (3.72) 

a d0 


w(x,0,z) - W°(0) (3.73) 

With the above restrictions a mathematical formulation equivalent to 
equations (3.32) to (3.36) becomes 


d0 


(3.74) 


D„ .4 o 
22 d w 


d0 


N. j 2 o 
0 d w 

a d0 


P " “ 


(3.75) 


N 


0 


A 1 2 e + A 22 


J_ dv° + + i. / i. dw° 

a d0 a 2 l a d0 


(3.76) 


Equations (3.74) to (3.76) constitute a sixth order system of ordinary 
differential equations for the dependent variables Ng, w°, and v°, 
with p and e as prescribed loads. They are solved subject to the 
six boundary conditions 


v°(±a) 


♦v 


(3.77) 


w°(± a ) - 0 


(3.78) 
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a 


dw . 
d9' (ta) 


(3.79) 


Equations (3.74) to (3.76), subject to boundary conditions (3.77) to 
(3.79), are solved by the same method of solution described in 
Section 3.3.5. The results of the solution are 


w°( 0) 



fi 2 . I P a f>, 6 | aa (cosh X9 - cosh Xa) 
W ; In \ sinh xe 


where 


(3.80) 




(3.81) 


and Ng and the pressure are related by 


A(N.) R 
6 


+ B(N 0 ) R 


+ c(N e ) 


(3.82) 


in which 
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pa - N 0 



A 



Xa \ + _1_ a 2 1_ a 3 

tanh Xa J 4 tanh Xa 4 2 

sinh Xa 
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Xa \ 
tanh Xa / 


(3.83) 
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(N e - a 12 '0 


22 
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v 

a 
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The parameter X In equation (3.81) is the magnitude of the nonzero 
roots (±X) of the characteristic equation for the homogeneous solution 
to w°(0), and represents the rate of exponential decay of the bending 
boundary layer. The root X of equation (3.81) is equivalent to the 
root Xj of equation (3.48) which is associated with circumferential 
shear and bending in the shear deformation theory. For the Klrchhoff- 
Love approximation there is no root equivalent to X 9 of equa- 
tion (3.48). A boundary decay length associated with circumferential 
bending may be calculated by substituting X from equation (3.81) into 
the definition given by equation (3.68). 

For geometrically linear strain-displacement relations equa- 
tions (3.74) to (3.76) of the Ki rchhof f-Love theory become 



D. . ,4 o 
22 d w 


d0 




(3.85) 


N 


e 






(3.86) 


These linear equations are a sixth order system of ordinary differential 
equations and are solved subject to the six boundary conditions given by 
equations (3.77) tc (3.79). Equation (3.85) may be easily Integrated 
and the solution derived is equivalent to that determined from beam 
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theory when the right hand side of equation (3.85) is eauated to an 
effective pressure Pj 


Pi 



(3.87) 


After integration of equation 
tions in equations (3.78) and 


w°(0) 




(3.8 5) 
(3.79) 



and application of boundary 
the solution for w° is 


condi- 


(3.88) 


The circumferential stress resultant Ng is still an unknown in equa- 
tion (3.88). The resultant Ng may be found by integrating equa- 
tion (3.86) subject to boundary conditions in equation (3.77) and using 
the fact that Ng is a constant from equation (3.84). The result of 
this operation is the following expression for Ng 


N 


e 



* 22 (t :)**# 


22 


1 + 


22 


4 5D 


2 4 
a a 


22 


2 4 

a a pa 


(3.89) 


The linear Kt rchhof f-Love solution for w° in equation (3.88) does not 
exhibit a bendirg boundary layer as does the nonlinear K1 rchhof f-Love 
solution given by equation (3.60) and the nonlinear shearing deformation 
solution given by equations (3.5G) to (3.60). ~Te linear and nonlinear 
Ki rchhof f-Love approximations are compared in Section 4.3. 
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3.4 Through-The-Thickness Stresses 

The transverse stresses (r , x , a ) are calculated from the 

9z xz zz 

elasticiry equilibrium equations (3.7) with the assumptions, consistent 
with the one-dimensional shell analysis, that all stresses are assumed 
to e independent ol the x-coordinate and the displacement w is 
i. it, erder.t of x as given by equations (3.26). With these assumptions 
equations (3.7) r iduce to 


>t . 

x t 

fte’ 



o 


(3.90) 


i .^00 


•> 0 


(3.91) 


. iee , i a f q &6 &w\ = 

r* dfc » n a 50\ a d0 / 


(3.92) 


T.xi cit »j».vessioni for the stress conponents a , 0 o , and r . 

x o x9 

in the thlcknes coo; di-ioc e follow from equations (3.10), (3.17), and 

(3.29), combined wii h , * C. Although o is not needed to deter- 

‘y ! x 

nune '.lie transverse sti-sse >n the elasticity equilibrium equations, it 
is give- •.•clow since it Is needed in subsequent failure analyses. The 
inplane stress components within a typical lamina are 
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Q,,e * 


Q, 
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i? r e 4 


i2 r eo 




(3.93) 


Qi 2 e 


g 22 e 0 + z 


' 2 ?* 00 


Q 26 F x0 


(3.94) 
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\e - Q 16 £ + Q 26 £ 0 + z 


Q 26 r ee + Q 66 r x6 


(3.95) 


The transverse shear stress x is determined by substituting x 
from equation (3.95) into the equilibrium equation (3.90) and integrat- 
ing in the z-direction. Similarly the transverse shear stress x^ z is 
obtained by substituting Oq from equation (3.94) into the second 
equilibrium equation (3.91) and integrating on z. The constants of 
integration resulting from these indefinite integrations are determined 
by requiring the transverse shear stresses to vanish at z = — t / 2 which 
is consistent with boundary conditions (3.8). Finally, the expressions 
for the transverse shear stress Xg z , the circumferential stress Oq in 
equation (3.94), and the deflection w°(0) in equation (3.52) are 
substituted into the third equilibrium equation (3.92) and indefinite 
integration on z is performed to determine a ? > The constant of 
integration is determined by setting o •= -p at z = -t/2. In the 
integration processes described it is important to recognize that the 
transformed reduced stiffnesses are piecewise constant functions of z. 
The results of this lengthy process are 


,, v , dT ft , dT 
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(3.96) 
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(3.101) 


In equations (3.99) to (3.101), z i designates the thickness coordinate 
to the top of the ith layer in the laminate; i.e., » (i - N/2) h, 



i = 1,2, .... N, where h is the layer thickness, and N is the number 

of layers in the laminate (t = Nh). If z is set equal to t./2 in 

equations (3.99) to (3.101), then q, (t/2) = A , q 0 (t/2) = B = 0, 

Umn mn’ ^2mn mn * 

and q ^ (t/2) = D , in which A , B , and D are the extension, 

Jmn mn ran mn mn ’ 

coupling, and bending stiffnesses, respectively, of classical lamination 

theory. The stresses x , x„ , and o given by equations (3.96) to 

xz dz z 

(3.98) vanish at the upper surface, as they should, because the shell 
equilibrium equations are integrated versions of the elasticity 
equilibrium equations. 

3.5 Strength of Materials Solution 

A simple strength of materials approach was used to perform a non- 
linear membrane analysis of long cylindrical panels under internal 
pressure. The strength of materials analysis gives good approximations 
for center deflections and circumferential membrane strains for 
sufficiently tbj.n panels. The circumferential membrane stress resultant 
can serve as an initial approximation for the shell stress resultant in 
the bending theory discussed in Section 3.3.5. Local effects such as 
the edge bending strain are not predicted from this analysis. An 
equivalent analysis for a curved plate is described in Reference 1. 

The basic assumption of this analysis is that the initially 
circular panel remains circular when internally pressurized. Under an 
internal pressure p the panel with initial radius a deforms to an 
arc of a circle with a radius R. This deformation process is shewn in 
Figure 6. Vertical equilibrium relates the circumferential stress 
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circular cylindrical membrane due to 


c>3 


resultant N p to the pressure p and deformed radius R by the 
relation 

N g - pR (3.102) 

The constitutive relationship depends on the assumed axial strain state 
as follows: 


For e ° = 0 
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e 22 e e 
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(3.103) 


These axial strain states are examined in Chapter 4. The geometric 
expression relating strain to deformation (change in radius) is 


o SR - aa 
E 0 aa 


(3.104) 


where the angles a and p are the semi-opening angles of the panel 
before and after deformation. Geometric compatibility requires 

R sin p = a sin a = C (3.105) 

where C denotes the half chord length shown in Figure 6. The 
equilibrium equation (3.102) is combined with an appropriate constitu- 
tive equations (3.103), the strain-displacement relation (3.104), and 
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geometric compatibility (3.105) to give an equation relating pressure to 
deformed radius. Using the constitutive relationship for the axial 
strain equal to zero, the equation is 


P 


R sin 


a sin 


C 

R 

C 

a 


1 . 


The solution may be obtained by solving this equation for the radius of 
the deformed panel R given the pressure p. Since this equation is 
transcendentally nonlinear in the unknown radius R, but linear in the 
pressure p, it is more expedient to solve the problem by choosing a 
value of R and calculating the corresponding pressure p. By succes- 
sive approximation the radius can be determined for any given pressure. 
With the radius determined, other quantities such as center deflection 
and edge rotation may be calculated from simple geometric considera- 
tions. The center deflection (w c ) and edge rotation ( * 0 - a) are 

w c = R[1 - cos(sin *(C/R))] - a[l - cos(a)] (3.106) 


y e = sin ^(C/R) - sin ^(C/a) 


(3.107) 


Chapter 4 


ANALYTICAL RESULTS 

In this chapter analytical results for the different analyses 
described in Chapter 3 are presented. In Section 4.1 two-dimensional 
finite element results are used to validate the one-dimensional approxi- 
mation solved analytically in Chapter 3. In Section 4.2 the geometri- 
cally nonlinear nature of the response of shallow cylindrical panels is 
illustrated by plotting center deflection and middle surface strain as a 
function of pressure for panels of differing stiffnesses. The circum- 
ferential surface strain distributions are given in Section 4.3 to 
illustrate the severity of the bending gradients at the edge, and to 
show the boundary layer decay lengths for the range of panels tested. 

In Section 4.4 the effects of a prescribed axial strain are discussed, 
and in Section 4.5 different boundary conditions are examined. The 
parametric studies in Sections 4.4 and 4.5 are useful for interpreting 
the experimental results, since ideal clamped boundary conditions were 
not achieved in the test fixture design. In Section 4.6 Kirchhoff-Love 
and shear deformation theories are compared to illustrate the importance 
of through-the-thickness shearing deformations in the bending boundary 
layer of the shell. Also the shear deformation theory for the one- 
dimensional analysis is used to assess the effect of the twist-curvature 
coupling ( D 25 t 0 ) of laminated composite shells. 
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4.1 Two-Dimensional Finite Element Results 

Geometrically nonlinear computer analyses with the finite element 
code STAGS (STructural Analysis of General Shells) (Ref. 18) were per- 
formed on two panels representative of those tested. The purposes of 
the STAGS analyses were to identify regions of the panels where a one- 
dimensional infinitely long cylindrical panel analysis wa6 applicable, 
and to validate the one-dimensional analytical solution. A 5-ply ortho- 
tropic panel and a 16-ply quasi-isotropic panel were chosen for the 
STAGS analyses because these panels represent extremes in thickness and 
lay-up of the test articles. Thus, they were expected to exhibit 
different response characteristics to the pressure loading. Both panels 
had a 60-inch radius, an 8-inch circumference, and were 20 inches long. 
Clamped boundary conditions were applied to all four edges. Since the 
STAGS computer code does not have shearing deformation capability, the 
clamped boundary conditions were modeled by constraining the three dis- 
placements along the edge and the slope normal to the edge to be zero. 

The 5-ply laminate had a layup of (±45/90/ +45] ^ where the sign con- 
vention for the fiber angle $ given in Figure 5 is used. A contour 
plot of the radial deflection is shown in Figure 7, and this figure re- 
veals that the deflection is spatially uniform in the x-direction over a 
center axial length of approximately 10 inches at a pressure of 1.3 psi. 
This central area is the region in which a one-dimensional analysis is 
appropriate. The length of 10 Inches for this central region of axial 
uniformity is slightly less than the 12-inch estimate for this geometry 
following the argument given by Douglas as repo^-d in Reference 1. 
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In Figure 8 the center deflection from the STAGS analysis is compared 
with center deflection from the one-dimensional shallow panel analysis 
as a function of pressure. There is excellent correlation between the 
two-dimensional and one-dimensional analyses. 

The circumferential surface strain distributions for the two- 
dimensional and one-dimensional analyses are shown in Figure 9 for the 
5-ply panel. Besides showing that excellent agreement is obtained for 
strains, Figure 9 illustrates an important advantage of the one- 
dimensional analysis. Strains may be determined at any circumferential 
point including the edge of the panel with the analytical solution. The 
strains in the STAGS analysis are determined at the location of the cen- 
troid of the finite element. Thus, interpolation is required to estimate 
strains at other circumferential locations in the finite element method. 
In particular, the strain gradient at the edge of the panel, which is 
important for failure analysis and in estimating through-the-thickness 
strains, may not be represented very accurately in the finite element 
solution since to compute it would require even higher derivatives of 
the basic displacement data. Since the majority of panels failed at the 
straight “dge, the one-dimensional analysis allows more insight into a 
critical aspect of the problem than the finite element solutions. 

Excellent agreement between the two-dimensional and one-ditnensional 
analyses was also obtained for the strains and displacements for a 
16-ply laminate. Because of the increased stiffness of this panel over 
the 5-ply panel, the portion of the axial length in which the one- 
dimensional analysis is valid is decreased substantially. This distance 
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Figure 0.- Normalized center deflection of a 5-ply r ±45/90] s panel from 
a two-dimensional nonlinear STAGS analysis and a one- 
dimensional nonlinear shallow panel analysis. 
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Figure 9.- STAGS and one-dimensional nonlinear shallow panel analysis results 
showing circumferential surface strain distribution along the 
circumferential center line (x-O.) of a 5-ply r±45/S0] panel with 
2.3 psi internal pressure. 8 
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Is illustrated in Figure 10. The distribution of axial surface strains 
indicates bending occurs for a distance of .6 to .7 X/L away from the 
curved edges. The axial strain is equal to zero in the center of the 
panel as was true for the 5-ply laminate. This axial membrane strain 
state resulted from requiring that the axial displacement along the 
straight edge of the panels be zero in the STAGS model. The axial 
membrane strain was treated cs a parameter in the one-dimensional 
analysis for infinitely long cylindrical panels. 

4.2 The Importance of Geometric Nonlinearity 

The stiffening effect caused by the geometric nonlinearity on the 
panel response with increasing pressure is sherm in Figures 11 and 12. 

In Figure 11 the pressure is plotted as a function of center deflection 
normalized by the shell thickness for a 5-ply [+45/ 75] s laminate having 
an initial radius of 60 inches and an 8-inch arc length. The results 
from two nonlinear analyses are shown; the strength of materials 
membrane solution of Section 3.5, and the Kirchhof f-Love shell theory 
(which includes bending) given by equations (3.80) to (3.83). In 
addition the linear Ki rchhof f-Love solutions given by equations (3.88) 
and (3.89) are also shown in Figure 11. As illustrated in this figure, 
there is little difference between the two nonlinear analyses which 
suggests that bending does not significantly affect the center deflec- 
tion. However, there are significant differences between the linear and 
nonlinear solutions shown in Figure 11. The nonlinear solution is 
considerably stlffer than the linear solution even at low pressures. 
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Figure 11.- A comparison of normalized center deflection response of a 5-ply 
[±45/90] predicted by three different theories. 
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For example, at a pressure of approximately 0.6 psi the linear theory 
predicts a center deflection equal to the panel thickness. The linear 
theory deflection is approximately 50 percent greater than the nonlinear 
theory deflection at 0.6 psi. 

Figure 12 also shows the stiffening effect of the geometric non- 
linearity on a plot of normalized pressure pa/A 22 versus circumfer- 
ential membrane strain for three different values of the semi- 

0 

opening angle a. These results were obtained from the nonlinear 
membrane analysis presented in Section 3.5 (bending is neglected). 
Shallow panels are completely characterized by the semi -opening angle 
o; smaller values of a imply shallower panels. Figure 12 illustrates 
that for a fixed value of the pressure a shallower panel has a smaller 
membrane strain. In addition, the curves in Figure 11 provide a good 
estimate of the middle surface circumferential strain at a given value 
of the pressure for the more complex nonlinear panel analyses, which 
include both bending and membrane actions, presented in Chapter 3. 
Consequently, the elementary membrane solution closely predicts the 
center deflection and circumferential strain, but by definition ignores 
the bending at the panel boundaries. 

The reason geometric nonlinearity is significant for a shallow 
internally pressurized panel can be found from the elementary membrane 
analysis in Section 3.5. For an increment in the arc length ds the 
following kinematic formulas can be derived from the equations in 
Section 3.5 for shallow panels: 
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dR 3 1 

ds “ ‘ 2 3 

a 

d_6 3_ i_ 

ds S 4a 2 
a 

dw - , 

c_ _3 _1_ 

ds 4 a ’ 

in which R is the radius, 0 is the rotation of the panel at its 
support, and is the center deflection. The above formulas are 

obtained with the chord length fixed. The initial radius and semi- 
opeui n t angle are denoted by a and a, respectively. These formulas 
show the decrease in radius, increase in edge rotation, and increase in 
center deflection are very sensitive to an increase in arc length fcr 
shallow panels. The increase in arc length may be a result of strain, 
or perhaps a slip at the edges in a test fixture. Thus there are sig- 
nificant changes in geometry in shallow panels for small changes in 
circumferential strain. Figure 13 shows the dramatic decrease in panel 
radius as the arc length increases for a panel with an initial radius 
a = 60 inches and initial arc length 2aa = 8 inches. Finally, since 
equilibrium of the circular membrane requires Nq = pR, it is clear that 
the decrease in radius ct a constant pressure in the nonlinear analysis 
results in a smaller circumferential stress resultant than obtained from 
a linear analysis. 

For laminates with considerable bending stiffness the nonlinear 
c.embrane effects becomes less important and bending resistance must be 
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taken into consideration. A solution to the linear shallow panel theory 
which allows bending and extension was examined. This solution was 
described in Section 3.3.6 (equations (3.87) to (3.89)). A nondimen- 
sional plot of membrane strain and center deflection as a function of 
nondimensional pressure is given in Figure 14 for values of the 
parameter \ defined in the figure. The 4- to 16-ply panels examined 
in this investigation have values of the X parameter which fall 
between values of 0 and 1. The lines representing an 8-ply and a 16-ply 
panel are shown in the figure. The 4-ply laminate would be represented 
by a line very close to the line \ = 0. The results given in Figure 14 
are accurate for cases where the center deflection of a panel are less 
than its wall thickness. For the 16-ply panel, the center deflection is 
greater than the wall thickness for pressures less than 10 psi. There- 
fore, linear results cannot provide accurate predictions for thin panels 
at pressures above 10 psi. 

The effect of panel thickness on the response is shown in 
Figures 13 and 16. Each panel has a common initial radius of 60 inches 
and a common middle surface arc length of 8 inches. The thicknesses of 
the 4-ply, 8-ply, and 16-ply panels are 0.0208, 0.0416, and 
0.0832 inch, respectively. The pressure versus center deflection in 
Figure 15 and the pressure versus membrane strain in Figure 16, are 
obtained from the nonlinear panel analysis of Section 3.3. These 
results show that the response is nonlinear even at moderate pressures 
fcr the thicker 16-ply panel. As would be anticipated, the thinner 


panels are more compliant 
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shallow panel theory results showing circumferential midd 
strain and nondlmensional center deflection response to 
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Figure 15.- Nonlinear shallow panel theory results showing center deflection 
response to pressure for a 4-, 8-, and 16-ply panel. 
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4.3 Circumferential Strain Distributions 

The circumferential surface strain distributions of an 8-ply 
laminate are shown in Figure 17 for the nonlinear membrane analysis 
(Sec. 3.5), and the linear and nonlinear shallow shell analyses 
(Sec. 3.3.6) with the Ki rchhof f-Love assumption. A great deal of varia- 
tion exists between these analyses. The nonlinear shallow cylindrical 
panel analysis exhibits an edge bending boundary layer visible in 
Figure 17 and represented mathematically by the exponential terms in the 
homogeneous solution given in equation (3.80). A boundary layer is not 
distinct for the linear solution given in Figure 17. The linear solu- 
tion given in equation (3.88) does not exhibit a bending boundary layer; 
i.e., the homogeneous solution does not have exponential functions in 
the solution. The strains predicted by the nonlinear shallow shell 
analysis are uniform away from the clamped edge. The bending strains 
increase substantially in a narrow region near the edge, which i6 
referred to as the edge bending boundary layer. The nonlinear membrane 
analysis approximates the center bending strain, but cannot predict any 
local edge bending. (The bending strains for the nonlinear membrane 
analysis were calculated from the change in curvature of the deformed 
and undeformed radii of the circular arcs shown in Figure 6.) Although 
the linear analysis overpredicto the maximum value of the tensile strain 
in this 8-ply example, it cannot in general be used for conservative 
design. In a 5-ply example, the edge surface strains determined from 
nonlinear analysis exceeded those predicted from the linear analysis on 



Nonlinear shallow shell 
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Figure 17.- Comparison between circumferential surface strain distributions 
generated by nonlinear membrane theory, and linear and nonlinear 
Kir chhof f -Love shallow shell theories. 
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both the concave and convex surfaces. The nonlinear and linear analyses 
agree more closely for thicker panels. 

Laminate bending stiffness plays an Important role in the cnaracter 
of the local bending boundary layer exhibited in the nonlinear shallow 
cylindrical panel analysis. The character of this boundary layer is 
represented in Figure 18 where normalized bending strain distributions 
are given for 4-, 8-, and 16-ply laminates at 20 psi pressure. The 
inside surface strains were normalized fo allow comparisons between the 
panels of three different thicknesses. The middle surface membrane 
strain was subtracted from the concave surface strain of each panel. 

The membrane strain is a spatial constant so the above subtraction 
yields the bending strain distribution. The distributions were 
normalized to go to the value of -1 at the center by dividing the bend- 
ing strain distribution by the value of center bending strain. 

Figure 18 shows the bending is uniform over most of the circumference 
for the 4-ply laminate with a boundary layer width of 0.75 Inches. The 
boundary layer width may be easily determined in the bottom portion of 
the figure where the vertical scale has been greatly magnified. From 
this enlargement it r av be seen that bending varies over the entire 
circumference of the 16-ply laminate. The edge bending 3train 
normalization factor is largest for the thinnest panel. The bending 
strain distribution for the convex surface would be Identical in shape 
but inverted with respect to the abscissa. 

The boundary dev_ay length of the Kirchof f-Love shallow shell theory 
is expressed mathematically by equation (3.68) where the value of X is 
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Figure 18.- Nonlinear shallow panel theory results showing normalized 

circumferential bending strein distribution on the concave surface 
of A-, 8-, and 16-ply panels. 
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given by equation (3.81). The boundary decay length is plotted in 
Figure 19 for 4-, 8- and 16-ply laminates as a function of pressure 
with e = 0.01 in equation (3.68). The boundary decay length is 
inversely proportional to the square root of the circumferential stress 
resultant and therefore decreases with pressure for all the laminates. 
Since the decay lengths are short and edge bending strains larger for 
the thinner panels, bending gradients are more severe for these panels 
than the thicker panels. 

4.4 Effect of an Applied Axial Strain 

The results of the previous sections are based on analyses which 
assume that the axial strain is zero. Tensile axial middle surface 
strains were measured experimentally and for this reason it is important 
to understand this effect. In this section the nonlinear membrane and 
bending analyses are used to examine three axial strain states. In the 
first case the axial stress resultant is set to zero. This state occurs 
when free axial contraction is allowed due to the Poisson affect and is 
the correct state for an arch as opposed to a long cylindrical shell. A 
negative axial strain results from the Poisson contraction. In the 
second case the axial strain is set to zero. This state is the one pre- 
dicted by the STAGS analysis in the central panel section away from the 
curved edges and is the correct state for a long cylindrical shell with 
axial displacements along the straight edge constrained. A tensile 
membrane stress resultant is present for the second case. In the third 
case a tensile membrane strain is applied. For purposes of this 
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comparison the value of this strain was chosen to be equal to the magni- 
tude of the negative strain which results from lateral contraction in 
the first case where the axial stress resultant is zero. Prescribed 
axial strains (e) and stress resultants (N x > were applied to an 8-ply 
quasi-isotropic [ ±45/ 9 0/ 0 ] laminate with a 60-ir.ch initial radius and 
8-inch arc length at a pressure of 20 psi. Numerical results are pre- 
sented in Table 2 for the membrane and bending solutions. Circumferen- 
tial membrane strain, circumferential stress resultant, and center 
deflection are given for the membrane solution in the table. The same 
physical quantities along with the edge bending strain and the middle 

surface interlaminar shear stress (x ) at the edge are given for the 

uz 

bending solution in the table. Both the circumferential membrane strain 
and the center deflection decrease as the axial strain increases for 
both the bending and membrane solution. Conversely, the circumferential 
stress resultant increases due to the Poisson effect. Comparing the 
membrane and bending theories .»e circumferential membrane strain and 
center deflection are approximately 9 percent and 2 percent larger, 
respectively, when calculated from membrane theory. For the bending 
theory the edge bending strain and shearing stress decrease as the axial 
strain increases. 

4.3 Circumferential Boundary Conditions 

Boundary conditions were examined to study their influence on both 
the panel center and edge responses. Also, the results were used to 
model closely the actual conditions on the panels as tested. This 
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Table 2. Effect of Axial Strain on Stress and Deflection Results 
Calculated with Membrane and Bending Theory (8- ply 
Specimen, p = 20 psi) 


Membrane Theory 



-.000623 0. .00196 630. 2.81 31.5 

0.0 206. .00181 648. 2.64 32.4 

.000623 412. .00167 667. 2.48 33.4 


Bending Theory 
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modeling of actual conditions was accomplished by measurement of 
boundary displacements. An additional reason for examination of 
boundary conditions is that the examination allows representation of 

actual fuselage condition? which are very difficult to simulate experi- 

o o 1 dw° 

mentally. Edge displacements, v and w , and edge rotation, — — — , 

a dB 

boundary conditions were examined and are discussed in individual 

1 dw° 

sections. These three conditions 'v°, w°, and — ) are given by 

v a dt) ' 

equations (3.77) - (3.79). The axial condition was satisfied by the 
assumptions that u° = 0 along the edge of the panel. 

4.5.1 Circumferential Displacement Condition 

The effect of circumferential displacements or slip is shown in 
Figure 20 for a 16-ply panel at 20 psi. In this figure the concave and 
convex edge strains are plotted as the circumferential slip is allowed 
to increase from 0 to 0.020 inches. The membrane strain is also plotted 
and goes through zero as the slip increases and becomes compressive for 
higher values of slip. The bending strains represented by one-half the 
difference between the concave and convex strains increased substan- 
tially in the process. These results are very informative in explaining 
events which occurred when the test panels slipped circumferentially in 
the clamped edge supports under internal pressurization. This circum- 
ferential slipping will be discussed further in Section 5.1. 

4.5.2 Radial Spring Boundary Condition 

In an actual fuselage the skin may deform radially at the stiff- 
ener. This deformation is difficult to duplicate experimentally without 
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testing a complete stiffened cylinder. However, it can be examined 
analytically by allowing the skin support to deform radially through the 
use of an extensional spring. To address the problem a special analvsis 
was conducted in which the boundary condition on w° (see equa- 
tion (3.78)) was replaced by the more general condition of a radial 
spring support with stiffness K. Shearing deformations were not 
included in this analysis. The solution procedure was similar to the 
procedure described in Section 3.3.5. The effect of radial spring 
stiffness on the membrane stress resultant for a 4- and 16-ply laminate 
is shown in Figure 21. As tne radial spring stiffness K goes to zero 
the membrane stress resultant approaches the complete cylinder hoop 
tension solution (hoop tension equals pressure times original radius). 
The deviation between the linear and nonlinear solution is larger for 
the 4-ply panel than the 16-ply panel. The 16-ply panel carries a large 
portion of the load through bending, whereas the 4-ply panel stretches 
more with the resulting decreased radius and carries the pressure 
through the nonlinear membrane action. For low radial spring stiffness 
these results indicate that geometrically nonlinear effects are less 
important. 

Circumferential surface strain distributions are shown in Figure 22 
at 30 psi pressure for one panel which has a radial spring stiffness K 
of 100 lb/in. and for one radially constrained panel. Results for these 
otherwise identical panels were obtained with the nonlinear shallow 
panel analysis. The difference in the membrane circumferential strain 
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Figure 21.- Effect of radial spring stiffness on the circumferential stress 



Concave surface 
Convex surface 
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Figure 22.- Circumferential surface strain distributions for 4-ply r±45l s 

panels which are rigidly and elastically restrained in the radial 
direction. 
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agrees with results shown earlier in Figure 21. The additional informa- 
tion determined from this plot is that severe bending gradients at the 
edge still exists for the finite radial spring stiffness. The largest 
strains still occur at the edgi, but the bending strains are decreased. 

4.5.3 Rotational Boundary Condition 

f 1 dw° \ 

The effects of edge circumferential rotation * - — -^-g-} are 
bound by considering clamped ( 4 / 0 * 0 ) and simply supported (Mg - 0) 
boundary conditions with shearing deformations suppressed. The solution 
for clamped boundary conditions is given by equation (3.80) and a solu- 
tion for simply supported boundary conditions is derived similarly. 
Radial, circumferential, and axial displacements (w°, v°, and u°) at 

the edge were set to zero to avoid the unnecessary introduction of addi- 
tional variables. 

The membrane strain response of an 8 -ply laminate is shown in 
Figure 23 for clamped and simply supported boundary conditions. The 
membrane strain is slightly higher for the simply supported condition. 
Although not shown, there is little difference between the center 
deflections for clamped and simply supported boundary conditions. The 
deflection for simply supported conditions are slightly greater than for 
the clamped condition. The difference between the membrane strain and 
central deflection for simply supported and clamped boundary conditions 
increases with bending stiffness. However, for the 4- to 16-ply panels 
teited, the rotational constraint did not have a large influence on the 
center deflection and circumferential middle surface strain. 


CO 



Figure 23.- Ciicumferential membrane strain response to internal pressur 
8-P- jr f±45] g panel with claaped and simply supported edges. 
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The edge bending strains were significantly influenced by the 
rotational edge support. Circumferential surface strain distributions 
are shown in Figure 24 for an 8-ply quasi-isotropic panel with simply 
supported and clamped boundary conditions. The bending strain is a 
maximum at the edge for the clamped condition and is zero for the 6imply 
supported condition. The predicted failure location for the clamped 
panel is at the edge and away from the edge for the simply supported 
panel. The rotation constraint therefore has a major effect on the 
response at the edges but '*nly a minor effect on the response at the 
center of the panel. 

4.6 lrfluence of Through-the-Thickness Shear Deformation 

Shear deformation is important for metal structures with a span-to- 
thickness ratio less than about 10. For materials which have a high 
Ei j /C 1 2 ratio, such as graphite-epoxy, shear deformation effects may be 
important for higher span-to-thickness ratios. The panels tested have 
span-to-thickness ratios greater than 80 but their graphite -epoxy mate- 
rial construction and the rapid variation of bending moment at the edge, 
make consideration of shear deformation important. The effects of 
through-the-thickness shearing deformations are approximated by allowing 
the normal to the middle surface to rotate independently of the circum- 
ferential and axial line elements in the deformed middle surface. Thus, 
the right angle between the normal and circumferential line elements is 

not preserved during deformation, i.e., shear strain y° i 0. 

z 6 



- Clamped 

- Simply supported 
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Figure 24.- Circumferential surface strain distributions for 8-ply [90/0/*45l 
p4J»el with simply supported and clamped boundary conditions. 
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Similarly, the right angle between the normal and the axial line ele- 
ments is not preserved during deformation, i.e., shear strain 

Y° ^ 0. For the long panel analysis developed in Section 3.3, the 
zx 

shear strain y° may also be nonzero since the normal can rotate 
'zx 

through an angle out of the z -0 plane even though the axial line 

element is restricted not to rotate. Hence v° = d> . However, 

t 0 occurs in the long panel analysis only if the twist coupling 
coef f icient °26 * U * In lamination theory D 26 j 4 0 reflects the 
material coupling between circumferential normal stress and the shear 


strain y on a parallel surface. The inplane shear strain y on a 
X 0 X 0 

parallel surface is determined, in part, by the change in the rotation 

of the normal out of the 9-z plane along the circumferential arc length, 

1 d *x 

or mathematically by — If f 0, as it is for the laminates in 

this study, circumferential bending is accompanied by a nonzero rotation 
1 34 x 

gradient — — — , which implies <\i f 0 , and hence a nonzero shear strain 
3 00 X 

Y results. 

'xz 

The distribution of the through-the-thickness shear strain y° 


predicted by the shear deformation solution of Section 3.3.5 is shown in 


Figure 25 for a 4-ply [+45] g laminate at 100 psi. The shear strain 
Y° is zero over most of the circumferential arc length in Figure 25 
except for a narrow region near the edge. In this edge region y° x 
changes very rapidly, attains a peak value very close to the edge, and 


then returns to zero at the edge. Since <l x is prescribed to vanish at 
the clamped edge, and y° = 4» v » Y° must vanish at the edge also. As 
discussed in the previous paragraph, if D 2 5 = 0 , then the shear strain 
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Figure 25. 
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Y zx would vanish over the entire arc length. Also the Ki rchhof f-Love 

shell theory presented in Section 3.3.6 would predict y zx = 0 over the 

o 

entire am length. The shear strain y zx i® zero because the rotation 
4 , x must vanish over the entire arc length to satisfy the hypotheses of 
the Ki rchhof f-Love theory for this problem, and not because D 26 = °* 

It is interesting to compare the predictions of the shear deforma- 
tion theory, designated SDTj, with the Kirchhof f-Love theory, desig- 
nated KLT. This is done in Table 3 where various response measures are 
computed from the two theories at pressures of 10 psi and 100 psi for 
4-, 8-, and 16-ply laminates. The table also presents results from a 
second shear deformation theory, designated SDT 2 , which was obtained 
from SDTj by setting °26 " °* Although setting D 2 ^ = 0 appears 
arbitrary, and is unnecessary in this analysis, setting D 2 ^ = 0 nust 
be done in more complicated two-dimensional laminated plate and shell 
analyses to make mathematical solutions tractable. For example Wang 
(Ref. 20) has shown the separation of variables technique does not work 
in solving the linear partial differential equations for laminated 
plates when and £> 2 ^ are nonzero. Separation of variables will 

work when and vanish in such problems (specially orthotropic 

plates, for instance). Since it is a relatively simple matter to set 

equal to zero in SDTj to obtain SDT 2 » this is done there to see 
the effect. 

Comparing the Kirchhof f-Love theory (KLT) to the first shear defor- 
mation theory (SDTj) for the laminates and pressures presented in 
Table 3 the following statements may be made: The center deflection 
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Table 3. Effect of Shearing Deformations on Panel Center and Edge 
Response Quantities 


Lamii'ate 


w 

c, 





and 

Pressure 

Aralysis 

v " t 

Uin. 

<r 9 ) e 

CD O 
N 

fl> 

Ct° ) 

xz < 

in. 

in. 


P 

P 


SDTj 

0.1303 

2097. 

-.5207 

-165.1 

34.69 

p*10ps-. 

sdt 2 

.1303 

2097. 

-.5001 

-142.2 

-46.79 


KLT 

.1303 

2097. 

-.5049 

-144.8 

-47.65 

r±«5] c 

SDT j 

.3924 

1066. 

-3.453 

-226.2 

47.55 

p=100psi 

sdt 9 

.3924 

1066. 

-3.208 

-194.8 

-64.14 


KLT 

.3921 

1065. 

-3.356 

-212.8 

-70.04 

[ -^+5 ] 2 S 

SDT] 

.08722 

1197. 

-.1366 

-63.5 

5.31 

p*10psi 

s^t 2 

.08722 

1197. 

-.1356 

-62.0 

-10.20 


KLT 

.08721 

1196. 

-.1364 

-62.7 

-10.31 

f±45] 2s 

SDT! 

.2874 

6449. 

-.9653 

-96.5 

8.40 

p=l OOpsi 

sdt 2 

.2874 

6449. 

-.9524 

-93.7 

-15.43 


KLT 

.2871 

6443. 

-.9803 

-99.0 

-16.30 

[iA5] As 

p=10psi 

SDT! 

SDT 0 

.0532 

.0532 

608. 

608. 

-.0385 

-.0385 

-30.7 

-30.5 

1.23 

-2.51 


KLT" 

.0532 

608. 

-.0387 

-30.7 

-2.53 

r*«) 4s 

SDTj 

.2033 

3699. 

-.2791 

-45.2 

1.88 

p=100psi 

sdt 2 

.2033 

3699. 

— . 2 "*8 3 

-44.9 

-3.69 


KLT 

.2030 

3693. 

-.2837 

-46.4 

-3.81 
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for SDTj is greater than w c for KLT. Circumferential middle 

surface strains e° are about the same for each theory. Differences 
9 

of less than 10 percent occur for the edge rotation gradient Tg 
between the two theories, but one theory does not consistently pre- 
dict larger magnitudes for Tg with respect to the other theory. 

Toe rotation gradient Tq is directly proportional to the circumfer- 
ential bending strain. Differences in the shear stress x° computed 

vZ 

on the middle surface at the edge are less than 20 percent for the two 

theories. Again one theory does not consistently predict a larger 

magnitude for x° with respect to the other theory. There are large 
9z 

discrepancies in the values of shear stress computed at the edge 

o 

of the middle surface. First, there is a sign difference in with 

KLT predicting negative values, and SDTj predicting positive values. 
Second, the magnitudes of from KLT are greater than the magni- 

tudes of t° from SDTi. These large discrepancies in x° are a 
xz i xz 

consequence of neglecting the rotation <1^ in Ki rchhof f-Love theory. 

The differences are largest for the thinnest panels. 

The effect of neglecting 1*2^ can be examined by comparing the two 
shear deformation theories SDTj and SDT 9 in Table 3. Both theories 
predict essentially the same center deflection and circumferential 

middle surface strain e°. Values of the edge rotation gradient T- 

e o 

differ between the two theories, SDTj predicts larger values than 
SDT 2 * The shear stress on the middle surface at the edge has a 

consistently smaller magnitude in SDT 2 with respect to SDTj. 

However, the shear stress x° at the edge of the middle surface has a 
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consistently larger magnitude in SDT 2 with respect to SDTj. What 

appears to be more important, though, is that SDTj predicts positive 

values of whereas SDT-, (like KLT) predicts negative values of 

o 

T xz* 

The comparisons between the three theories given above do not 
establish the correctness of any one over the others. If an exact non- 
linear elasticity solution were available to compare to, then a con- 
clusive statement on which theory is best may be made. No exact 
elasticity solution is available to compare with the three theories, 
however. Since the first shear deformation theory is less restrictive 
in its assumptions relative to the other two theories, it is assumed to 
be the best of ti e three. In Chapter 5 the theory used to compare to 
experiment is SDTj which will be referred to as the shear deformation 
theory. In the following section of this chapter further comparisons of 
the three theories are made to illustrate differences in the predicted 
distributions of o 2 , t 0z , and t xz through-the-thickness of the 
laminate. 

4.7 Through-the-Thickness Stresses 

In this section the through-the-thickness stress components t , 

0Z 

t , and 0 are examined in more detail. The interlaminar stresses at 
xz z 

the panel middle surface (z*0) and along the circumferential arc from 
the center to the edge of a panel are shown in Figure 26 for [±45] 
symmetric laminates of 4 and 16 plies. These stresses are calculated at 
an applied internal pressure of 100 psi and are normalized by this 
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Figure 26.- Circumferential distribution of normalized interlaminar stresses at 
middle surface of 4- and 16-ply T±45] symmetric panels with 100 psi 
internal pressure. Results calculated with SDT. theory. 
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applied pressure. For the 4-ply panel the interlaminar stresses are 

nearly zero over most of the arc length except for a small distance near 

the edge. For the 16-ply panel the interlaminar stresses are uniform 

for about ore-half of the distance to the edge. The shear stresses 

and x° are zero at the center with the shear stress 
9z xz 9z 

o 

increasing monotonically toward the edge. The slope of the shear 

o 

stress distribution changes sign near the edge. The normal stress a z 
is compressive except in a small region a short distance away from the 
panel's edge where it is tensile. The maximum values of all inter- 
laminar stresses occur at the edges of the panel, and the maximum values 
are larger for the thinner panels than for the thicker panels. Also, at 

the panel edges, the magnitude of the shear stress is greater than 

0z 

the shear stress which, in turn, is greater than the normal stress 

o°. The severe interlaminar stress gradients shown in Figure 26 occur 
z 

near the panel edges as a consequence of the severe bending strain 
gradient shown in Figure 18. The interlaminar stresses are proportional 
to the bending strain gradient as shown in equations (3.96) to (3.98). 

The results shown in Figure 26 are obtained from the first shear 
deformation theory (SDTj) which includes the effects of the twist 
coupling coefficient ^26* figures 27 and 28 these interlaminar 
stresses are also computed from Ki rchhof f-Love theory (KLT) and the 
second shear deformation theory (SDT 2 ), and plotted for comparison 
with SDTj. In Figure 27 the results plotted are for a 4-ply [±45] 
symmetric laminate at 100 psi, and in Figure 28 the results are plotted 
for a 16-ply (±45] symmetric laminate at 100 psi. The distributions are 
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Figure 27.- Circumferential distribution of normalized interlaminar 

stresses calculated at the middle surface of a 4-ply r±45l 
panel with 100 psi internal pressure. Results calculated 
with KLT, SDTj, and SDT 0 theories. 
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Figure 28.- Circumferential distribution of normalized interlaminar 
stresses calculated at the middle surface of a 16-ply 
[±45]^ s panel with 100 psi internal pressure. Results 
calculated with KLT, SDT^ , and SDT 9 theories. 
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only plotted near the edges ( 0/ a * 0.90-1.0). For both laminates the 

largest differences occur at the edge with no differences occurring away 

from the edge. The most pronounced difference between the three 

theories occurs for the interlaminar stress x° . In the mathematical 

xz 

solution the roots ±X^ given by equation (3.48) are responsible for 

these differences. These roots are associated with the rotation 4> x 

and the twist coupling coefficient ^26* differences in the 

stresses t° and o° between the theories are greatest at the edge 
9z z 

but small in magnitude. Through-t he-thickness shear and twist coupling 
effects extend over a greater distance from the edge for the 16-ply 
panel. 

The boundary layer decay lengths from shear deformation theory 
(SUTj) are associated with roots Xj and equation (3.48), and 

are given in Table 4 for 4-, 8-, and 16-ply (±43] symmetric laminates. 
These decay lengths are calculated with equation (3.68) after evaluating 
the roots with equation (3.48). The value of the root an order 

of magnitude larger than the root Xj for the pressures examined. The 
root Xj is numerically close to the root X (eq. (3.81)) of the 
Kirchhof f-Love theory and both Xj and X increase with pressure and 
their corresponding boundary decay length decreases with pressure. For 
the Internal pressures loads which caused small strains in the panels, 
the boundary decay length associated with X 2 is much shorter than the 
decay length of Xj. In addition, the boundary decay length associated 
with X 2 remains approximately constant with pressure. Both boundary 
decay lengths increase with increases in laminate thickness. 
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Table 4. Boundary Decay Lengths and ^ 4- » and l&~ply 

[±45°] Laminates with Internal Pressures of 10 and 100 psi 

ns 




* 

* 

Laminate 

p. psi 

in. 

Cp * ^ • 

r±45] s 

10 

0.6223 

0.0586 

100 

0.2730 

0.0582 

' ±45 >2s 

10 

1.638 

0.1372 

100 

0.722 

0.1369 

f ± 4 S U s 

10 

100 

4.582 

1.883 

0.2534 

0.2833 


ill 


The distributions of stress components o , x Q , and z through 
the thickness are shown In Figure 29 for the 1+45] symmetric laminates 
of 4 and 16 piles. These results are calculated at the clamped edge 
using SDTj. The through-the-thi ckness shear stresses are symmetric 
about the middle surface and the normal stress is nonsymmetric. The 
shear stress x. is the largest in magnitude relative to x and a 

W XZ Z 

and attains a maximum magnitude at the middle surface. The shear stress 

x^ is a maximum at the clamped edge of the middle surface. The 

maximum values of the stresses x and o occur at locations other 

xz z 

than the middle surface. 


The influence of the twist coupling coefficient on the 

through-the-thickness distributions of stress a z , Xg z , and x xz is 
shown in Figures 30 and 31 for the same 4-ply and 16-ply [ + s5] symmetric 
laminates. The twist coupling coefficient has the largest influence on 
the shearing deformation stress x . The shear deformation theory 
(SDTj) predicts the largest value of the stress x^ in Figure 30(a) 
for the 4-ply panel whereas tne Ki r chhof f-Love theory predicts the 
largest value in Figure 31(a) for the 16-ply panel. Therefore, as 


pointed out previously, the Ki rchhof f-Love theory does not, in general, 


give conservative solutions. 

In this paragraph through-the-thickness stresses of two quasi- 
lsotropic 8-ply symmetric laminates and one [ ±4 5 ] 8-ply symmetric 
laminate are con^ared. The through-the thickness distribution of o z , 
x^, and x xz for three 8-ply laminates of different laminate 
stacking sequences are shewn in Figure 32. These distributions -*re 
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Figure 29.- Through-the-thickness distribution of normalized 

interlaminar stresses at clamped edge of 4- and 16-ply r ±45 ] 
symmetric panels with 100 psi internal pressure. Results 
calculated with SDT theory. 
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Figure 30.- Through-the-thickness distribution of normalized interlaminar 
stresses at clamped edge of a 4-ply [±45] s panel with 100 psi 
internal pressure. Results calculated with KLT, SDTj, and SDT 2 
theories. 




Figure 31.- Through- the- thickneFs distribution of normalized 

interlaminar stresses at clamped edge of a 16-ply [±45] 
panel with 100 psi internal pressure. Results 
calculated with KLT, 5DT, , and SDT. theories. 
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Figure 31.- Concluded. 
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Figure 32.- Through-the-thickness distribution of normalized 

interlaminar stresses at clamped edges of three different 8- 
ply panels with 100 psi internal pressure. Results 
calculated with SDT theory. 





1 17 


calculated at the clamped edge using the shear deformation theory 

(SDTj). The laminate stacking sequence has a large influence on all 

through-the-thickness stresses. The largest change in slope of the 

shear stresses x and Xq_ occurs at the interface of ±45° plies, 
xz fte 

The largest slope changes for the normal stress (o z ) distributions occur 
at the interface of the 90° ply and 0° ply on the concave side of the 
middle surface. 



Chapter 5 


COMPARISON OF ANALYSIS AND EXPERIMENT 

To correlate the analysis with the experiment, it was found 
necessary to measure the initial shape of the panel after fitting it 
into the test fixture, and to measure the slipping of the panel out of 
the test fixture with increasing pressure. Using the measured initial 
geometry and the measured pressure-slip characteristics in the analyti- 
cal model, the comparisons presented in this chapter establish that the 
analysis can predict reasonably well the experimental response data. 

With confidence established in the analytical model of the experiment, 
it is reasonable to use the stresses and strains predicted in the 
analytical model to analyze failure, which is the subject of Chapter 6. 

Results of thirteen test specimens are given in Table 1. Most of 
the comparisons between analysis and experiment are carried-out on 
specimens G2, G4, G7 , G9, and G10 listed in Table 1, because only these 
panels were surveyed prior to pressurization to determine their actual 
shape after fitting into the test fixture. The other panels were tested 
before it was determined that the initial shape deviated substantially 
from a nominal 60-inch radius cylindrical panel. Thus, comparison of 
analysis and experiment is presented only for the panels In which the 
actual shape was determined and used in the analysis. In particular, 
the comparisons are presented for the 4-ply (G2), 8-ply (G7), and 16-ply 
(G10) panels. Nc comparisons are presented for the 5-ply (G4) and 
16-ply (G5) panels, since these results are similar to those presented 
for specimens G2, G7 , and G10. 
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5.1 Measurements Required for the Analysis 

5.1.1 Panel Stiffnesses 

Some of the elastic stiffnesses utilized in the analysis were 
verified by actual testing. Tensile coupons were cut from flat panels 
of the same material and cured under the same conditions as the 5-ply 
specimens G3 and G4 and the 8-ply quasi-isotropic specimens G5, G6, and 
G7. These panels were tested in tension to failure. The specimen 
elastic stiffnesses were measured and compared to their values computed 
from classical lamination theory using the lamina properties given in 
Section 2.1 for T300-5208 and the average lamina thickness measured from 
the coupons. The experimentally determined extensional stiffnesses of 
the 5-ply specimens (G3 and G4) were 97,900 lb/in. and 164,200 lb/in. 

The theoretically predicted stiffnesses are 103,500 lb/in. and 

167.000 lb/in., respectively. The experimentally determined extensional 
stiffness of the quasi-isotropic 8-ply specimens (G5, G6, and G7) was 

325.000 lb/in. and the theoretical stiffness was 321,500 lb/in. The 
experimental and theoretical stiffnesses agree very well and suggest 
that the use of lamination theory with the elastic properties for 
T300-5208 given in Section 2.1 is sufficient for predicting the elastic 
stiffnesses of the specimens. The computed stiffnesses used in the 
analysis for panels G2, G4, G7, G9, and G10 are given in Table 5. The 
average lamina thicknesses used in computing the stiffnesses are also 
given in Table 5 for each panel. 


Table 5. Laminate Thicknesses and Elastic Stiffnesses Used in Correlation of Analysis 
and Experiment for Specimens G2, G4, G7, G9, and CIO. All Laminates Are 
Symmetric and the Coupling Stiffnesses =0. 
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5.1.2 Initial Geometry 

Shape surveys of unpressurized specimens clamped in the test 
fixture were conducted to determine the actual initial geometry of the 
specimens. The instrumentation and procedures used to conduct these 
initial geometry surveys are described in Section 2.2. The results of 
typical measurements to determine the shapeb of circumferential (x - 0) 
and axial (6 = 0) lines passing through the center of specimen G7 are 
shown in Figures 33 and 34, respectively. A circular arc with a radius 
of 39.9 in. has the same rise and chord as the measured shape and is 
shown in Figure 33 for comparison. Points on the measured curve fall 
inside of the circular arc and the radius of the circular arc is 33 per- 
cent smaller than the nominal 60 inch panel radius. The radius of the 
circular arc was used in the analysis to approximate the effects of the 
initial panel geometry. The same procedure was used to determine radii 
of 33.4, 43.1, 53.4, and 50.3 inches for panels G2, G4, 09, and G10, 
respectively. The shape of the measured axial line through the panel 
center shown in Figure 34 differs from the nominal straight line genera- 
tor of a cylindrical panel indicating that the initial shape of the 
specimen is also curved in the axial direction. The curvature in the 
axial direction was neglected in the analyses. These differences in 
shape indicate that the initial geometry of the panel is distorted from 
the nominal circular cylindrical panel shape. 
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Figure 33.- Measured vertical shape along the circumferential center 
line (x«0.) of specimen G7 before pressurization. 
Appropriate circular shape used in analysis shown for 
comparison . 
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Figure 34.- Measured vertical shape along the axial center line (0*0.) of 
specimen G7 before pressurization. 
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5.1.3 Pa lie 1 Slip 

Applied internal pressure loading causes circumferential and axial 
tensile forces to be generated in a specimen. These tensile forces can 
c '.use the specimen to pull away or slip from the test fixture clamping 
mechanism and this slipping affects the response of the specimen. 

In-plane circumferential displacements at the specimen edges were 
measured with direct-current differential transformers a6 the applied 
internal pressure was increased to determine the magnitude of displace- 
ments or slippage which occurred during testing. Results of typical 
ln-plane circumferential displacement measurements at the midlength of a 
straight edge of specimen G7 are shown in Figure 35. The pressure- 
displacement curve in Figure 35 is approximately bilinear up 'o failure. 
Circumferential displacement values from these pressure-displacement 
curves were used as prescribed values for v (eq. (3.37)) in the one- 
dimensional panel analysis to approximate the effects of ln-plane cir- 
cumferential edge displacements on panel response. Rotations at the 
clamped edges also occur, but were difficult to measure for inclusion 
into the analysis. 

The two-dimensional STAGS analysis results shown in Figure 10 indi- 
cate that the axial strains are zero at the center of an ideal panel 
with perfectly clamped edges. However, in the experiment back-to-back 
axial strain gage6 located at the center of the specimens indicate that 
tensile axial membrane strains are nonzero. The axial strain curves for 
specimen G7 are shown in Figure 36. The 4-, 5-, and 8-ply specimens 
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Figure 35.- Measured circumferential edge displacements u. ed in analysis 
of panel G7. 
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tested in this study had nonzero axial membrane strains at the specimen 
center. The stiffer lb-ply specimens had both nonzero axial membrane 
and bending strains at the specimen center. The nonzero axial membrane 
strains are a result of the two-dimensional effects associated with 
pillowing and the inability of the test fixture clamping mechanism to 
prevent the specimen from slipping axially at the straight edges of the 
panels. In the one-dimensional analysis, the assumed axial boundary 
displacements of equation (3.25) are directly responsible for the 
uniform axial strain state of equation (3.29). Although the axiaJ 
boundary slip was never measured, axial slip probably occurred simul- 
taneously with the measured circumferential slip. This assumption is 
substantiated by comparing Figures 35 and 36 which exhibit a change in 
slope at essentially the same pressure (27.6 psi) for specimen G7. 
Therefore, the two-dimensional axial strain effects, caused by slipping, 
are modeled as an applied axial strain e in equation (3.65) of the 
one-dimensional analysis. 

5.2 Comparison of Analysis and Experiment 
5.2.1 Center Deflection Versus Pressure 

The effects of increasing internal pressure p on the nondiroen- 
sionalized radial deflections at the center w c are shown in Figure 37 
for the 4-, 8-, and 16-ply panels (specimens G2, G7, and G10, respec- 
tively). The radial deflection is nondimensionalized by the panel 
thickness. The solid curves in the figure represent results from the 
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Figure 37.- Experimental and analytical radial deflection response at 
centers of specimens G2, GV and G10. 
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one-dimensional panel analysis and the symbols represenc experimental 
data. The filled symbols represent specimen failure. The effects of 
the initial measured radius (e.g. , Fig. 33), the measured in-plane 
circumferential edge displacements (e.g., Fig. 35) and measured axial 
strains at the panel center (e.g.. Fig. 36) are included in the analysis 
of each panel. Analytical results for the 8-ply panel (G7) that do not 
include the effects of measured initial radius, in-plane circumferential 
edge displacements and axial strain at panel center are shown in 
Figure 37 by the dashed curve. Comparing the results of the analysis 
for the 8-ply panel with the experimental results indicates that the 
effects of initial geometry and panel slip must be included in the 
analysis for accurate correlation with experimental results. The non- 
linear character of the global pressure-displacement responses in 
Figure 37 is indicated by the large radial deflections relative to the 
panel thicknesses and the increases in the slopes of the pressure- 
displacement curves, or stiffening, as internal pressure increases. 

These analytical results correlate well with the corresponding experi- 
mental results up to failure. 

5.2.2 Membrane Strain Versus Pressure 

Circumferential membrane strains at the centers of the typical 4-, 
8-, and 16-ply panels (specimens G2, G7, and G10, respectively) as a 
function of internal pressure are shown in Figure 38. The solid curves 
represent the analytical results and the symbols represent the experi- 
mental results. The filled symbols represent specimen failure. 
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Figure 38.- Experimental and analytical circumferential membrane strain 
responses at center of specimens G2, G7 and G10. 
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These results indicate that the circumferential membrane strains at the 
panel center are larger at every pressure for the thinner panels than 
for the thicker panels, and that the analytical results for the two 
thinner panels correlate well with the experimental results up to 
failure. Although the analytical results for the 16-ply panel correlate 
reasonably well with the experimental results for the lower values of 
applied internal pressure, the analytical and experimental results do 
not agree as well at the higher values of the internal pressure. The 
back-to-back axial strain gages at the center of the 16-ply panel 
indicate that negligible bending strains exist at the lower pressures, 
but the strain gage data show that both bending and membrane strains are 
significant at the higher pressures and increase as the pressure is 
increased. For the thinner panels the bending strains were negligible 
up to failure. Apparently, axial bending effects at the center are more 
important for the stiffer 16-ply panels than for the thinner panels and 
become more pronounced as the applied pressure is increased. 

5.2.3 Circumferential Strain Distributions at a Fixed Pressure 

Circumferential surface strain distributions at the panel midTength 
(x = 0, Fig. 1) are shown in Figure 39 for the typical 4-, 8-, and 
16-ply panels (specimens G2, G7, and G10, respectively) discussed in the 
previous section. The distributions are shown for only half the panel 
due to symmetry and are for an applied internal pressure of 50 psi. 
Results from the one-dimensional panel analysis are represented by the 
solid curves and the experimental results are indicated by the symbols. 



= 50 psi 
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Figure 39.- Experimental and analytical circumferential surface strain 
distributions for an internal pressure of 50 psi. 
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Strain gage data for the concave surfaces are represented by the 
circular symbols and data for the convex surfaces are represented by the 
square symbols. These results indicate •’hat the circumferential 
membrane strain found by averaging the two surface strains is nearly 
constant along the circumferential arc length for each panel as was 
predicted by the one-dimensional panel analysis. The surface strains 
are uniform near the center of the panels, but significant bending 
strains exist near the panel edges. These bending strains are maximum 
at the edges. The magnitudes of the bending strains are larger for the 
thinner panels than for the thicker panels, and the bending strain 
gradients are confined to a smaller edge region for the thinner panels. 
For example, the bending strain of the 4-ply panel G2 are distributed 
over approximately 10 percent of the circumferential arc length (see 
Fig. 39(a)), but are distributed over about 80 percent of the arc length 
of the 16-ply panel G10 (see Fig. 39(c)). The strains predicted by the 
analysis agree with the experimental strains and indicate the severity 
of the bending strain gradient near the panel edges. The magnitudes of 
the strains measured near the edges are smaller than those predicted by 
the analysis. These differences may be due to slight rotations of the 
specimens and clamping bars at the panel edges which relieve the bending 
strain. Also, differences may occur due to the difficulty of measuring 
strains in a region of severe bending gradient. 
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5.3 Prediction of Through-the-Thickness Stress Components for 

Specimens G2 and G10 

The through-the-thickness stresses c z , , and Tq z were cal- 

culated from the analyses of specimens G2 (4-ply) and G10 (16-ply) and 
the results are presented in Figures 40 and 41 for an applied internal 
pressure of 50 psi. The stresses shown in these figures are normalized 
by this applied pressure. The distribution of interlaminar stresses 
along the circumferential arc and at the panel middle surface for the 
4-ply specimen G2 is shown in Figure 40(a) and those for the 16-ply 
specimen G10 are shown in Figure 40(b). The normal stress a z is 
compressive except for a small tensile region near the panel edge. At 
the panel edges the magnitude of the shear stress is greater than 

the shear stress r xz , which is greater than the normal stress o z * The 
maximum values of all interlaminar stresses are at the panel edges and 
the maximum values are larger for the thinner panels than for the 
thicker panels. The bending strain distribution in Figure 39 and the 
equilibrium equations (3.90) to (3.92) indicate that the thinner panels 
with severe bending gradients near the panel edge (e.g., see Fig. 39(a)) 
also have severe interlaminar stress gradients near the panel edge 
(e.g., see Fig. 40(a)). The distributions of through-the-thickness 
stresses at the panel edge for the 4-ply panel are shown in Figure 41(a) 
and those for the 16-ply panel are shown in Figure 41(b). These results 
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Figure 40.- Interlaminar stress distributions along the 
z-0, and at 50 psi for panels G2 and G10 on 
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Figure 41.- Through-the-thickness distribution of normalized interlaminar 
stresses at the clamped edge (0-ci) and for an internal pressu 
50 psi for panels G2 and G10. 
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indicate that the shearing deformation stresses are symmetric about the 
middle surface and the transverse normal stress is nonsymmet ri c. The 
normal stress is equal to the applied internal pressure on the concave 
surface of the panel and equal to zero on the convex surface. The shear 
stresses are both equal to zero on the panel external surfaces. The 
fluctuations in the shear stress x xz are related to changes in ply 
orientation from +45 to -45, and the maximum value of shear stress x 

0Z 

occurs at the panel middle surface. 



Chapter 6 


FAILURE 

6.1 Observations of the Failed Specimens 

6.1.1 Examinations of the External Surfaces 

All specimens except the aluminum panel A2 were tested to failure 
to study their failure characteristics, and to determine their ultimate 
pressure. The applied internal pressure was slowly increased until the 
specimens ruptured along a straight or curved edge and the applied 
pressure could not be maintained. The failures initiated at the bound- 
aries of the panels and not in the interior. Most failures occurred 
along the straight edges as shown in Figure 42. Photographs of the 
concave and convex surfaces of each failed panel are included in 
Appendix A. 

The edge failure locations for each panel shown in the photographs 
in Appendix A are summarized in Table 6. The straight edges are identi- 
fied as SI and S2 and the curved edges as Cl and C2. Edge SI is closest 
to the strain gages located along the circumferential center line, and 
edge Cl is closest to the strain gages located along the axial center 
line. The edge failures have visible damage on both the concave and 
convex surfaces. When the failure is characterized by a complete break 
of all lamina6 through the thickness, the edge is marked with an aster- 
isk superscript in Table 6. In the other cases the laminate failure 
does not completely separate the interior test portion from the clamped 
portion of the laminate. For example, the concave surface plies of 
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Figure 42.- Typical failed specimens. 



Table 6. Eoge Failure Description 



Panel 

Stacking 

Sequence 

Edge Failure 

Concave Ply Failure 
Only 


A1 

— 

S2* 

— 



A2 

— 

Not failed 

— 



Cl 

U45] s 

Cl* S2* C2 

SI 



G2 

l*«5] 8 

Cl S2* C2 

— 



03 

r ±45/90] 

8 

Cl*-S2*-C2* 

— 



G4 

[ ±45/90] 

s 

SI* C2* S2* 

— 



G5 

[90/145/0] 

s 

SI* 

C2 



G6 

[90/0/145] 

s 

SI* 

Cl 

S2 

C2 

G7 

[90/0/145] s 

S2* 

Cl 

SI 

C2 

G8 


Cl*-S2* 

— 



G9 

[145/145/90-/0,] 

L Is 

S2 

— 



G10 

[145/145/90, /0,] 
l is 

SI 

— 



G 1 1 

[90 2 /l45/i45/0,] s 

SI* 

— 





^An asterisk superscript on edge indicates all lamina have failed 
through- the- thickness 

2 

A dash between edge designators indicates the failure is continuous 
between edges 
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panels G9 and G10 delaminated and failed at the bolt holes as shown in 
Figures All and A12 in Appendix A. In Table 6 a dash between two edges 
indicates that the edge failure is continuous between these two adjacent 
edges. Also in the fourth column of Table 6, the edges in which the 
concave surface plies failed are identified. These edges did not fail 
catastrophically but the concave surface ply failures may have contrib- 
uted to the catastrophic failure of the other edges. De laminations 
occurred at the failed edge(s) of all the specimens. 

The 4- and 5-ply specimens Gl, G2, G3, and G4 were composed pri- 
marily of 45° plies and failed with extensive surface fiber delamination 
which extends from the edge to the center as shown in Figures A3, A4, 
and A5 in Appendix A. These panels had more damage around the curved 
edge chan the thicker panels. In general the thicker panels failed 
exclusively along the straight edge. The concave surface plies of 
specimens G6 . nd G7 have failed at the edge as indicated in Table 6. 
These failed plies are delaminated from the panel almost completely 
along the edge. Some delaminations extend as much as a distance of 
2 inches normal to the edge. The 16-ply specimens G9 and G10 had 
failures which extended to the bolt holes on the concave surface but not 
completely through the thickness of the specimen at the edge. Specimens 
G9 and G10 had the ±45° plies located closer to the outer surfaces of 
the laminate rather than near the middle. The 16-ply panel Gil which 
had ±45° plies at the middle surface failed completely through the 


thickness at the edge 
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6.1.2 Panel Dissection 

Typically the failures occurred along one edge of the panel leaving 
the other edge intact. The failed edge was often too severely damaged 
to determine where failure initiated. The intact edge opposite the 
failed edge of the panel experiences an identical stress state as the 
failed edge prior to ultimate failure. These stress states are iden- 
tical because of the symmetry in the pressure load and boundary condi- 
tions. Therefore, the intact edge should give a good indication of the 
local failure state immediately before failure. The results described 
in Table <* are from a surface examination of the specimens and describe 
the edge failure locations, surface damage states and whether the crack 
is continuous through the thickness or not. However, it was impossible 
to accurately determine what was happening through the thickness of the 
laminate. To accomplish this examination, the panels were cut along the 
circumferential direction (at x = 0) and the cross section was 
polisheu. Only the circumferential cross section will be examined since 
most failures occurred along the straight edge where the one-dimensional 
analysis was valid. The thinner panels had a greater tendency to fail 
near the curved edges. The circumferential cross section was examined 
under a 10 power hand lens and photographs were taken to determine the 
extent of any local damage at the intact edge opposite the failed edge. 

Photographs of the circumferential cross section of the intact and 
failed edges of the 16-ply specimen C9 are shown in Figure 43(a). 

Matrix cracks in the center four 0° plies are visible at both edges and 
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(a) 16-ply specimen G9 intact and failed edge cross sections 



(b) 8-ply specimen G7 intact edge cross section 


Figure 43.- Cross sections of intact and failed specimen edges. 
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may have been caused by the high transverse shear stress and the 

circumferential normal stress Og shown in Figure 41(b). A photograph 
of the cross section of the intact edge of 8-ply quasi-isotropic speci- 
men G7 is shown in Figure 43(b). Transverse matrix cracks and delamina- 
tions are visible, and the two outer plies on the concave side have 
failed in tension and have delaminated. 

The results of these cross-sectional examinations are summarized in 
Table 7. Damage at both the intact and failed edges are described in 
Table 7. The type of damage and the through-the-thickness ply locations 
are identified. Delaminations (D) and fiber breakages (F) of specific 
plies are identified in the table. The damage designation (D cr F) is 
followed by numbers in the table vhich either identify the adjacent 
plies which have delaminated or identify the ply which has fiber 
failures. The plies are consecutively numbered from the concave surface 
to the convex surface. Transverse matrix cracks were visible in most of 
the specimen cross sections and are not included in Table 7. It is very 
interesting that panels with significant failure along one edge had no 
damage along the intact edge. 

6.2 Interpretation of Strain Gage Data 

In addition to the photomicrographs the response of strain gages 
located at the clamped edge of the panels was also very important in 
understanding the failure mechanisms. The strain response of back-to- 
back circumferential strain gages near the straight edge of specimen G7 
is shown in Figure 44. The responses of the convex and concave gages 
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Table 7. Observable Damage In Machined Circumferential Edge 
Cross Section 




DAMAGE 


SPECI- 

STACKING 

INTACT 


FAILED 


MEN 

SEQUENCE 

EDGE 


EDGE 

COMMENTS 

G1 


D 3 l-2 b 


Severe 

D under tab 

G2 

[±451 

D 1-2 


Severe 

D 1-2 under tab 


s 

D 3-4 



D 3-4 in test 






section 

G3 

[±45/90l s 

D 1-2 


Severe 

D in test section 

G4 

[±45/5tf] g 

Severe 


Severe 

Both edges failed 

G5 

[90/145/0] 

None 


D 1-2 

D 1/8" from edge 





D 2-3 

D between all 





D 3-4 

plies at failed 





D 7-8 

edge 

G6 

[ 90/0/ ±45] 

F a l and 

2 

D 2-3 



s 

D 2-3 
D 3-4 


D 3-4 


G7 

[90/0/145] 

F 1 and 

1 

D 2-3 

D mainly in test 



D 2-3 


D 3-4 

section 



D 3-4 


D 5-6 


G8 

( !45 bs 

D 5-6 


D 3-4 

Both edges intact 
at cross section 
examined 
D under tab 

G9 

[±45/145/90-/0,] 

None 


D 6-7 

D 6-7 7/8" under 


Z Z S 



D 10-11 

tab and 2" in test 





D 13-14 

section 






Severe damage at 
failed edge 

G10 

[±45/145/90-/0,] 

D 5-6 


D 6-7 

Most damage on 



D 6-7 


D 10-11 

intact edge under 



D 2-3 



tab;F 1 and 2 under 
tabs 

Gil 

[90-/±45/±45/0-] 
l l s 

None 


D 2-3 

Most D in test 




D 5-6 
D 6-7 

section 



D = delamination; F = fiber breakage 

Numbers refer to lamina location. Lamina numbered consecutively 
from concave to convex surface. A dash indicates delamination 
occurs between laminas (i.e. 1-2 indicates delamination between 
the first and second ply) 
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8-PLY STRAIN RESPONSE 



• wv • vv/y v/ • • vvv • vw / 

CIRCUMFERENTIAL STRAIN, e 0 


Figure 44.- Circumferential strain response of back-to-back gages 
located at clamped edge (x-0., 6 - of specimen G7. 
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are smooth until a small discontinuity occurs at approximately 61 psi. 
This discontinuity is associated with an abrupt circumferential slip at 
the boundary shown by the displacement response data in Figure 35. The 
effect of circumferential slip on the edge strain response was shown in 
Figure 20. The slip caused an increase in bending strain and a decrease 
in middle surface strain. The decrease in the middle surface strain is 
offset by the increase in the bending strain at the concave surface. 

The opposite effect occurs at the convex surface and the slip shows up 
as a discontinuity in the gage response. 

Figure 44 also shows a larger discontinuity at an internal pressure 
of 82.7 psi where both the convex and concave strains abruptly change. 
The concave gage which was indicating a large tensile strain before the 
discontinuity registers zero immediately after the event. The convex 
gage which was increasing in compression registers tension immediately 
after the event. As the pressure increases the concave gage changes 
only slightly and continues to register nearly zero and the tensile 
strain measured by the convex gage increases. This behavior is 
explained by a local failure of the concave surface plies at the clamped 
straight edge. This concave surface failure is substantiated by the 
photomicrographs of specimen G7 shown in Figure 43. The 90° concave 
surface ply in Figure 43 has failed by a tensile fiber failure at the 
edge. The 0° inside ply next to the 90° ply has also failed circum- 
ferentially in a tensile matrix failure mode at the edge. Delamination 
of both the 0° ai.d 90° plies is apparent from the edge inward. This 
delamination propagates under the concave surface strain gage which is 



located approximately one-eighth of an inch away from the edge. Tne 
strain gage on the concave surface remains attached to the two delami- 
nated plies and registers a strain of zero reflecting the loss of load 
carrying capability. With the two plies on the concave side failed, the 
remaining laminate must now support the tensile stress resultant, and the 
convex strain gage registers an increasing tensile strain. The re- 
sponses of the edge strain gages for panels G2, G4, G7, G9, and G10 are 
given in Appendix B. Local edge failures, such as occur at 82.7 psi for 
specimen G7, have a large detectable influence on the strain gage 
responses shown in Appendix B. 

6.3 First Major Damage Event and Ultimate Failure 

The discontinuities in the edge strain gage response are believed 
to be caused by either a fiber fracture or delamination. Many local 
failures were audible during the test and their effects on panel 
response were detected by the edge strain gages. As will be shown in 
Section 6.4, the analysis predicts transverse matrix cracking to occur 
at lower pressures than fiber fracture. Transverse matrix cracking is 
classified as minor damage since it did not cause a discontinuity or 
abrupt slope change in the measured strain response. The strains at the 
first major damage event, and at ultimate pressure were determined from 
the edge strain response given in Appendix B and are presented in 


Tables 8 and 9 
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Measured circumferential surface strains near the clamped straight 
edge of the panel are given in Table 8 and axial strains near the 
clamped curved edge of the panel are given in Table 9. These strains 
are recorded by edge gages immediately before and after the first major 
damage event and at the ultimate pressure. The axial gages are located 
approximately at the center of the curved edge, and the circumferential 
gages are located at the center of the straight edge. The gages were 
either located one-eighth of an inch away from the edge or as close as 
possible to the edge which is indicated by 0 + in Tables 8 and 9. For 
example, distances of 0.020 in. and 0.030 in. were measured between the 
strain gages and the edge for specimens G9 and G10, respectively. The 
strains could not be recorded to failure for all panels because in some 
instances the strain gages failed due to the high strains. Also, for 
some of the panels in the table, the strains after the first major 
damage event are not reported since the first major damage event 
occurred simultaneously with rupture at the ultimate pressure. 

The first major damage event pressure and the ultimate pressure are 
plotted in Figure 45 for the specimens tested in this investigation. 

Also shown on the figure is the maximum pressure applied to the thicker 
aluminum specimen (specimen A2) which was not tested to failure. All 
specimens failed at pressures well above 20 psi which is greater than 
the nominal pressure used for ground test verification of pressurized 
fuselage structure. The aluminum specimens yielded at the specimen 
edges to form plastic hinges that reduced the bending strains at the 
edges. The lighter weight graphite-epoxy specimens failed at lower 
pressures than the aluminum specimens with the same nominal thicknesses. 
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ORIGINAL PAGfc IS 
OF POOR QUALITY 



Figure 45.- First major local failure pressure and ultimate pressure of 
all specimens tested. 
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Ply orientation and laminate thickness appear to have a strong 
influence on the ultimate pressure capability of the graphite-epoxy 
specimens. As shown in Figure 45, the 16-ply quasi-isotropic panels 
(G9, G10, and Gil) did not fail at twice the ultimate pressure of the 
8-ply quasi-isotropic panels (G5, Gb, and G7) and the 8-ply ±45°-angle- 
ply panel (G8) did not fail at twice the ultimate pressure of the 4-ply 
±45°-angle-ply panels ( G 1 and G2). Although the strains, bending 
gradients and interlaminar stresses at the panel edges are more severe 
for the ti. inner panels (see Figs. 39-41), the thinner panels are more 
structurally efficient on a weight basis with respect to their ultimate 
pressure capability than the thicker panels. The 8-ply quasi-isotropic 
panels (G5, G6, and G7) failed at lower pressures than the 8-ply +45°- 
angle-ply panel (G8) and the 5-ply 1+45/90]^ panels (G3 and G4) failed 
at lower pressures than the 4-ply +45°-angle-ply panels (Gi and G2). 
Graphite-epoxy specimens Gl, G3, G5, G9, and Gil ruptured with a sudden 
loss of pressure and with no indication from the response of the edge 
strain gages of progressive local failures or damage occurring at 
pressure below the ultimate pressure. The other graphite-epoxy speci- 
mens had one or more major local failures (e.g. , fiber fracture or 
delamination) occur at pressures below the ultimate pressure. 

6.4 Evaluation of Failure Criteria 

The stresses computed from the one-dimensional panel analysis were 
used to determine if the first major damage event could be predicted by 
some commonly used failure criteria. Stresses at the straight edges of 
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panels G2, G4, G7 , G9, and G10 were determined at the pressures asso- 
ciated with the major damage events given in Table 8. The analyses 
utilized the measured radius, circumferential slip, and axial strain for 
each specimen to model the response as accurately as possible. The 
strains and stresses determined frum the analyses were substituted into 
phenomenological and individual mode failure criteria. Six different 
individual modes of failure and two phenomenological failure criteria 
were examined. The individual failure modes included tensile fiber, 
compressive fiber, tensile matrix, compressive matrix, tensile inter- 
laminar, and compressive interlaminar failure modes (Ref. 21) and are 
shown in Appendix C. The material strength properties used in the 
failure criteria are given in Table 10. The phenomenological failure 
criteria include a two-dimensional Tsai-Wu failure criterion which is 
based on the in-plane stresses (o xyI °09» an£ * T x0^’ anc ^ t ^* e more 
general three-dimensional Tsai-Wu functional which includes the in-plane 
stresses and the through-the-thickness stresses (t xz > Tq z , and a„ z ). 
The Tsai-Wu functional utilizing the full three-dimensional stress state 
is given in Appendix D. The two-dimensional Tsai-Wu failure criterion 
found in many composite material text books is a specialized case of the 
three-dimensional criterion. These failure criteria were examined 
across the circumference and through the thickness of the panels. The 
maximum value (dimensionless) of each criterion occurs at the clamped 
edge but at different locations through the thickness. These values are 
given in Table 11 at the first major damage event pressures. Values 
greater than unity indicate failure in each criterion. The value of the 
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Table 10. Material Strength Properties of T300/5208 Used in the 
Failure Criteria 



„ 1 
Property 

Symbol 

Strength, 

ksi 


Axial Tensile Strength 

X t 

218.0 

Axial Compressive Strength 

X 

c 

-218.0 

Transverse In-plane lensile Strength 

Y t 

6.1 

Transverse In-plane Compressive Strength 

Y 

r 

-21.4 

Transverse Thickness Tensile Strength 

z 

t 

6.1 

Transverse Thickness Compressive Strength 

z 

c 

-21.4 

Shear Strength in x-y Plane 

S 

9.8 

XV 


Shear Strength in x-z Plane 

S ' 
xz 

9.8 

Shear Strength in y-z Plane 

S 

yz 

6.0 


Coordinate x refers to the axial or fiber direction 
Coordinate y refers to the transverse in-plane direction 
Coordinate z refers to the transverse thickness direction 
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Tade 11. Maximum Values of Failure Criteria at First Damage 

Pressure . All Maximum Failure Criteria Values Occur 
at the Clamped Edge.^ 


Failure Criterion 


Specimen 

Pressure, 

psi 

2-D Tensor 
Polynomial 

3-D Tensor 
Polynomial 

Tensile 

Fiber 

Compressive 

Fiber 

G2 

61.93 

22.41 (I) 

22.40 (I) 

1.29 (I) 

0.69 ( 65 " 

G4 

86.69 

27.45 (I) 

27.44 (I) 

1.44 (I) 

0.93 (0) 

G7 

82.55 

7.12 (1-2*) 

7.56 (1-2*) 

1.35 (I) 

0.66 (0) 

G9 

75.94 

5.64 (I) 

5.63 (I) 

0.64 (4-5*) 

0.42 (0) 

G10 

99.69 

3.75 (2-3*) 

3.78 (2-3*) 

1.11 (I) 

0.83 (0) 


Failure Criterion 


Specimen Pressure, Tensile Compressive Tensile Compressive 

psi Matrix Matrix Interlaminar Interlaminar 


G2 

61.93 

44.04 (I) 

10.99 (I) 

2.45 (M) 

2.29 (M) 

G4 

86.69 

55.09 (I) 

13.70 (I) 

3.67 (2-3*) 

3.45 (3-4*) 

G7 

82.55 

15.27 (1-2*) 

1.69 (1-2*) 

0.82 (M) 

0.76 (M) 

G9 

75.94 

8.97 (I) 

2.32 (I) 

0.45 (M) 

0.40 (9-10*) 

G10 

99.69 

5.41 (2-3*) 

1.43 (2-3*) 

0.53 (M) 

0.46 (9-10*) 


The through-the-thickness failure prediction location are identified by the 
following notation: 

I - Failure predicted at the inside or concave surface fiber 
0 - Failure predicted at the outside or convex surface fiber 
M - Failure predicted at the middle surface 

1-2* - Failure predicted at the interface of plies 1 and 2, and in ply 2 
marked with an asterisk. Plies are numbered consecutively 
beginning with the concave surface ply 
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tensile matrix lailure mode is much greater than unity which indicates 
that tensile matrix cracking was predicted at a pressure lower than the 
pressure of the first major damage event. The two-dimensional and 
three-dimensional Tsai-Wu functions are also very large due to the 
influence of a similar transverse matrix term in the criteria. No 
abrupt changes in the panels response nor substantial water leaking was 
observed due to the occurrence of transverse matrix cracking. The 
maximum value of the three-dimensional Tsai-Wu function occurs at the 
same location as the two-dimensional Tsai-Wu function and the values of 
each do not significantly differ. The influence of the through"the- 
thickness stresses which are a maximum at or near the middle surface is 
not significant on the magnitude or location through the thickness of 
the maximum Tsai-Wu value. The values of the failure criterion do 
change significantly at the midpiane of the laminate where the through- 
the-thickness stresses are largest and the bending stresses vanish. For 
most of the panels the tensile fiber mode has values closer to unity 
than any of the other criteria. For four of the panels the analysis 
predicts tensile fiber mode of failure values greater than unity. The 
analysis did not model all phenomena which were observed in the experi- 
ment. Transverse matrix cracking and edge rotation from incomplete 
clamping affect the experimental response but are not modeled analyti- 
cally. The tensile fiber mode is analogous to the maximum-stress 
failure criterion applied in the tensile fiber-direction and appears to 
correlate reasonably well witn the pressures corresponding to the first 
major damage event. 



Chapter 7 


CONCLUDING REMARKS 

An experimental and analytical investigation was conducted to study 
the nonlinear response and failure characteristics of internally pres- 
surized 4- to 16-ply-thick graphite-epoxy cylindrical panels with 
clamped edges. The panels were selected to simulate the skin between 
two frames and two stringers of a typical transport fuselage structure. 
The stiffness of each panel tested, as measured by the slope of the 
response curve relating the radial deflection of the panel center to the 
applied internal pressure, increased with increasing internal pressure 
which is characteristic of a geometrically nonlinear response. Clamping 
the panel edges caused local bending and interlaminar stress gradients 
near the panel i-dges and these gradients were found to be more severe 
for the thinner panels. The radial deflections of the panels were uni- 
form in the axial direction over a substantial portion of the central 
region of the panels. 

A one-dimensional cylindrical panel analysis based on nonlinear 
shallow shell theory was derived and used to determine the nonlinear 
response of the central region of the panels. The formulation included 
through-the-thi ckness shear deformations which introduce the twist- 
coupling coefficient into the analysis. Twist-coupling coeffi- 

cient effects were determined to be important for response quantities in 
the region of severe bending gradients near the clamped straight edge. 
Neglecting through-the-thickness shear effects in the analysis gives 
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edge stresses which are higher than those predicted with through-the- 
Lhickness shear for some laminates, but lower for other laminates. This 
nonconservative nature limits the use of the Kirchhof f-Love shell 
theory. When the measured initial radius and the effects of measured 
in-plane circumferential displacements at the panel edges and measured 
axial membrane strain at the panel center are included, the one- 
dimensional panel analysis accurately predicts the nonlinear response 
of panels away from the curved edges. The analytical results correlate 
well with the experimental results up to the first major damage event 
including the severe local bending gradients at the panel straight 
edges. 

Different boundary conditions were examined to study the nature of 
the geometric nonlinearity, to approximate more closely "real life" 
conditions, and to evaluate experimental boundary condition anomalies. 
The important geometric nonlinearity in this problem results from the 
product of the membrane hoop tension and the slope of the radial deflec- 
tion. The nonlinearity disappears when the edge is allowed free radial 
expansion. In this case the radial deflection slope goes to zero since 
every point on the circumferential curve deforms uniformly in the radial 
direction. Another boundary condition evaluated was the rotation. The 
circumferential and radial boundary displacements were set to zero and 
the limiting cases of clamped and simply supported edge rotation were 
examined. The differences between clamped and simply supported boundary 
conditions in the panel responses at the panel centers were very small 
for thin panels. Rotations of the panel edges allowed by incomplete 
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clamping and individual ply failures at the edge did not have a large 
influence on the center response of the panels. The fact that the panel 
center response was relatively insensitive to the extremes in rotational 
constraint imposed by clamped and simply supported boundary conditions 
explains why such good correlation existed between the test and clamped 
analysis results for pressures up to the ultimate pressure. The local 
bending strains near the edge were vastly different for the two dif- 
ferent boundary conditons. The simply supported panels had no bending 
strains at the edge. The clamped panels had rapid variations in the 
bending strain near the edge with the maximum occurring at the edge. 

The graphite-epoxy panels failed along the panel edges where the 
local bending and interlaminar stress gradients occur. Some graphite- 
epoxy panels ruptured with a sudden loss of pressure and with no indica- 
tion of local failures or damage occurring at pressure below the ulti- 
mate pressure. The other graphite-epoxy panels had one or more major 
local failures at pressures below the ultimate pressures. The failures 
appear to be caused by tensile lamina failures at the panel edges. 
Transverse matrix cracking and delaminations also occurred in regions 
with local bending and interlaminar stress gradients. Aluminum panels 
tested for comparison yielded and formed plastic hinges at the panel 
edges before rupture. A nonlinear analysis is required to predict 
accurately the stresses in regions with severe local bending and inter- 
laminar stress gradients. These accurate local stress predictions are 
necessary in order to predict the onset of failure in brittle graphite- 
epoxy laminates. All graphite-epoxy panels failed at pressures well 
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above the 2C psi nominal pressure used for ground test verification of 
pressurized fuselage structure. The experimental results indicate that 
ply orientation and laminate thickness have a strong influence on the 
failure characteristics and ultimate pressures of graphite-epoxy curved 
panels. 

In conclusion, the present study has identified and resolved 
important aspects in the response to internal pressure of an aircraft 
fuselage skin. The response was determined from experiment to be non- 
linear. A one-dimensional geometrically nonlinear analysis was 
developed and correlated well with the measured panel response. Both 
testing and analysis identified the skin region adjacent to the stiffen- 
ers as the critical area of the pressurized panel. A bending gradient 
exist in the skin near the stiffeners and initiated failure of all 
panels tested. The edge gradient response of the aluminum panels was 
different than the response of the composite panels. The aluminum 
relieved the edge bending gradient by yielding, whereas, the brittle 
graphite-epoxy developed transverse matrix cracks, fiber breakage, and 
delaminations at the edge. Although all panels carried pressures much 
greater than the design ultimate pressure, the long term durability 
issues of pressurized composite skins must be examined. 
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Figure AI.- Aluminum Specimen A1 
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Figure Am.- Graphite-epoxy Specimen C2; ft45 
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Figure A5.- Graphite-epoxy Specimen G3, [±^5/§0] 
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Figure A7.- Craphite-epoxv Specimen C5; [90/±45/0] g 




Concave Surface 


Figure A8.- Graphite-epoxy Specimen G 6; [90/0/iAS] 
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Figure A9.- Graphite- epoxy Specimen G 7 ; f90/0/±^5] 
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Figure A10.- Graphite-epoxy Specimen G8; f±45] 
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Figure All.- Graphite-epoxv Specimen G9; [±45/±45/90_/0 o ] 
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Figure A 1 2 . - Graphite-epoxy Specimen 
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Figure Bl.- Circumferential strain response of back-to-back gages 

iocated 0.030 inches from the straight edge of specimen 
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Figure B2.- Axial strain response of back-to-back gages located 0.040 
inches from the curved edge of specimen G2. 









Figure B4.- Axial strain response of back-to-back gages located 0.050 
inches from the curved edge of specimen G4. 










Figure B7.- Circumferential strain response of back-to-back gages 

located 0.020 inches from the straight edge of specimen 
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Figure B8.- Axial strain response of back-to-back gages located 0.125 
inches from the curved edge of specimen G9. 
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Figure BIO.- Axial strain response of back-to-back gages located 
0.125 inches from the curved edge of specimen G10. 
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Appendix C. Failure Mode Criteria 

The following failure criteria are used to predict different 
failure inodes. Values equal to or greater than one indicate failure in 
all the criteria. These failuie criteria were extracted from Reference 
[ 21 ]. 
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Appendix D. Three-Dimensional Tsai-Wu Failure Criterion 


The Tsai-Wu failure criterion is a special case of the tensor 
polynominal failure criterion where only linear and quadratic terms 
are included. The Tsai-Wu failure criterion is commonly used for 
two-dimensional stress states but can be expanded to three-dimensional 
stress states. For three-dimensional stress states the Tsai-Wu 
failure criterion has the following form 


F 1 C 11 + F 2°/2 + F 3°33 + F ll°ll + 2F 12°11°22 + 2F 13°1 1°33 


+ F-^o _ + -F__a 79 o - + F__a _ + 4 


22 22 
+ F 66 T 1 2 


23 22 j 3 33 33 4 44 l 23 ' ‘55 l 31 

1. (D.l) 


In the above expression the coefficients F^ , F^ and F^ which 
correspond to the linear x 7 3 , T 3 j and ^^2 terms are zero * These must 
be zero since the shear strength is independent of the loading 
direction. The coefficients in the failure criterion are determined 
in terms of the strength properties given in Table 10 by considering 
one-dimensional loadings. Consider the case of a tensile loading in 
the fiber- or one-direction such that 

Ojj 4 0 and o 22 = o 33 * t 23 * t 31 - t 12 * 0. (D.2) 

At failure Ojj * X t and the Tsai-Wu failure criterion becomes 


F i x t + F n x t * 1 


(D.3) 


Next compressive loading in the fiber direction is considered where 


at failure o 


11 


X and the Tsai-Wu failure criterion becomes 
c 
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F.X + F. ,X‘ = 1 
1 c 11c 


(D.4) 


The strength coefficients Fj and Fjj are determined by solving Eqs. 
(D.3) and (D.4) simultaneously. 


F - + J_ 

X, X 
t c 


(D.5) 


(D. 6) 


Similarly it may be easily shown that 


F 0 J + l_ 

Y Y 

t c 


(D.7) 


(D.8) 


from consideration of tensile and compressive loadings in the transverse 
or two-dire<~tion • Furthermore two additional coefficients may be 
determined from loadings in the thickness direction. 


1 + 1 


(D.9) 
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(D. 10) 


By considering the three possible shear loadings separately the following 
coefficients may be determined 
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(D. 11) 



194 


and 



(S ) 

xy 


2 


The coefficients Fp, F ?3 and may only be determined from biaxial 
stress states. Such tests are very difficult to conduct and these 
coefficients were assumed to be equal to the following values 
F 12 = F 03 - F 13 * -0.58 x 10" 10 


This value was found to give good correlation with off-axis tests for 
boron/epoxy by Pipes and Cole [22]. when used for F, 9 in the two-dimensional 
Tsai-Wu criteria. 
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